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ABSTRACT 

A  geochemical  baseline  is  being  determined  for  sur- 
ficial  materials  in  the  vicinity  of  oil  shale  tract  C-a  in 
the  Piceance  Basin  of  Northwestern  Colorado.   Two  hundred 
fifty-three  composite  samples  of  A-horizon  soils,  big  sage 
(Artemisia  tridentata  tridentata  and  subspecies  wyomingensis) , 
rice  grass  (Oryzopsis  hymenoides) ,  and  wheat  grass  (Agropyron 
smithii)  were  collected  on  a  grid  eight  by  six  miles. 

The  elements  studied  are  B,  Mo,  Zn,  Li,  As,  and  Hg. 
Concentrations  range  from  crustal  averages  for  Zn,  Hg,  Mo, 
and  Li  to  three  to  six  times  the  average  for  B  and  As.   Due 
to  alkaline  conditions  these  elements  may  pose  special 
environmental  problems  during  the  development  of  oil  shale 

resources . 

Analysis  of  variance  techniques  were  used  to  determine 
those  elements  that  display  regional  variations.   Hypothesis 
tests  were  employed  to  determine  the  significance  between  the 
means  of  components  found  on  the  Parachute  Creek  member  of 
the  Green  River  Formation  and  Uinta  Formation.   Most  com- 
ponents are  substantially  enriched  in  the  surficial  materials 
found  on  the  Parachute  Creek  member.   Regional  trends  are 
displayed  in  the  surficial  materials  of  the  Uinta  Formation 
for  several  components,  thus  regional  variations  are  not 
caused  exclusively  by  lithologic  changes.   There  are  no 
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discernible  differences  between  the  subspecies  of  sage  for 
the  elements  analyzed.   Composite  sampling  reduces  low 
scale  variance  for  elements  that  show  regional  trends  but 
local  variance  is  dominant  for  elements  whose  distribution 
is  fairly  homogeneous. 
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INTRODUCTION 

The  renewed  interest  in  oil  shale  resources  of  Colorado 
has  generated  concern  in  the  potential  environmental  impact 
of  surface  and  in- situ  retorts  in  Colorado,  Wyoming,  and 
Utah.   Numerous  studies  are  in  progress  to  determine  the 
regional  geochemistry  of  these  areas.   The  chemistry  of 
soils  (Klusman  and  Ringrose,  1976b),  streams  (McNeal  et  al, 
1976),  and  plants  is  studied  to  establish  baselines  for 
future  reference  in  monitoring  the  effects  of  the  develop- 
ment and  environmental  impact  of  the  oil  shale  industry. 

Oil  shales  are  fine  grained  rocks  which  contain  sub- 
stantial amounts  of  organic  material  that  can  be  refined 
into  fuels.   The  organic  material  is  separated  into  two 
fractions.   The  soluble  bitumen  fraction  makes  up  20  percent 
of  the  organic  material  and  the  remainder  exists  as  a  complex 
insoluble  material  termed  kerogen.   The  shale  oil  is 
obtained  by  "retorting"  the  pulverized  shale.   In  the  sur- 
face operations  the  rock  is  subjected  to  temperatures  of 
500-550°C.   This  breaks  the  chemical  bonds  holding  the 
organic  material  to  the  rock  matrix  and  vaporizes  the 
organic  material.   It  is  then  condensed  into  a  crude  oil 
that  can  be  refined  into  final  products  (Yen  and  Chilingarian, 
1976). 
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Surface  retorting  requires  the  mining  and  storage  of 
large  amounts  of  rock  and  waste  material.   The  oil  shale 
available  on  tract  C-a  is  best  obtained  by  surface  mining 
because  of  the  proximity  to  the  surface  of  the  richest 
deposits.   Southeast  from  tract  C-a,  oil  shale  tract  C-b  is 
also  being  evaluated.   On  oil  shale  tract  C-b  the  shale  is 
too  deep  for  economical  surface  mining  so  underground  or 
in-situ  methods  of  shale  oil  extraction  are  being  examined. 
Both  retorting  methods  produce  reducing  and  alkaline  con- 
ditions which  mobilize  many  harmful  trace  elements. 

Northwestern  Colorado  includes  parts  of  the  Southern 
and  Middle  Rocky  Mountains,  Wyoming  Basin,  and  Colorado 
Plateau  provinces  as  defined  by  Feneman  (1931) .   The  area  is 
drained  by  the  Yampa,  White,  Colorado,  and  Gunnison  Rivers, 
all  of  which  are  westward  flowing.   The  Piceance  Basin  which 
contains  an  estimated  1,200  billion  barrels  of  oil 
equivalent  (Murray  and  Haun,  1973)  is  the  major  structural 
feature  of  the  Colorado  Plateau  Province.   During  the 
deposition  of  the  sediments  in  the  Piceance  Basin  the 
climate  was  semi-tropical  with  flora  similar  to  that  found 
now  along  the  Gulf  of  Mexico  (Bradley,  1963) .   The  present 
climate  is  semi-arid  and  has  changed  the  chemical  conditions 
in  the  sediments  making  many  of  the  trace  elements  more 
easily  mobilized. 
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Purpose  and  Scope 

This  thesis  investigates  the  spatial  distribution  of 
trace  elements  in  surficial  materials  in  the  vicinity  of  oil 
shale  lease  tract  C-a  in  the  Piceance  Basin  of  northwestern 
Colorado  (Figure  1).   The  primary  objectives  of  this  study 
are:   1)  to  determine  the  average  values  and  ranges  of 
trace  elements  in  the  various  media,  and  2)  to  determine  the 
sources  of  the  natural  variations  in  the  geochemical  environ- 
ment.  Other  goals  of  the  survey  are  to  evaluate  the 
applicability  of  the  grid  sampling  design  and  develop  geo- 
chemical maps  for  those  components  which  display  statis- 
tically valid  variations. 
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Figure  1.  Simplified  geological  map  of  the  Piceance  Basin 
{Ringrose,1976  ) 
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GEOLOGIC  AND  GEOCHEMICAL  SETTING 

The  Piceance  Basin  is  a  large  northwest  trending 
structural  downwarp  bounded  on  the  north  by  the  White  River, 
on  the  east  by  the  White  River  Uplift  and  the  Grand  Hog- 
back, on  the  south  by  the  West  Elk  Mountains  and  Uncompahgre 
Plateau,  and  on  the  west  by  the  Douglas  Creek  Arch. 

The  Piceance  Basin  was  an  area  of  deposition  during 
the  Eocene  epoch.   This  fact  is  shown  by  the  thickening  of 
sediments  towards  the  axis  of  the  basin  and  thinning  towards 
the  margins  of  the  basin.   The  time  of  formation  of  the  basin 
was  probably  during  early  Eocene  time  because  of  the  presence 
of  basal  Wasatch  sediments  thickening  towards  the  center  of 
the  basin  (Donnell,  1961).   At  the  end  of  the  Cretaceous 
and  during  Tertiary  time,  the  Laramide  and  the  post  Laramide 
orogeny  formed  the  basis  for  the  structural  features  seen  in 
the  Piceance  Basin  today.   After  the  deposition  of  the 
extensive  oil  shales  of  the  Green  River  Formation,  many 
basalt  flows  were  extruded  on  the  broad  flat  area  of  the 
Piceance  Basin.   Eocene  and  post  Eocene  structural  deforma- 
tion were  the  most  intensive  episodes  of  folding  affecting 
the  region.   Further  uplift  in  middle  Pliocene  time  caused 
the  5000  feet  of  downcutting  expressed  by  steep  cliffs  in 
the  area  today.   More  gentle  uplift  continued  throughout 
the  Quaternary  period  creating  the  more  subdued  landscape 
in  the  central  portion  of  the  basin  (Murray  and  Haun,  1974). 
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The  most  important  stratigraphic  units  are  shown  in 
Figure  3.   The  Wasatch  Formation  is  a  varicolored  claystone 
and  shale  displaying  shades  of  gray,  red,  purple,  and  green. 
It  also  has  sections  of  massive  to  crossbedded,  mostly  len- 
ticular sandstone.   The  age  is  considered  to  be  at  or  near 
the  Paleocene-Eocene  boundary. 

The  Green  River  Formation  overlies  the  Wasatch  Formation 
and  contains  the  richest  deposits  of  oil  shale.   The  Green 
River  Formation  is  divided  into  three  members.   The  Douglas 
Creek,  Garden  Gulch,  and  Parachute  Creek  members  consist  of 
massive  to  platy  dolomitic  marlstone.   The  Mahogany  zone  of 
rich  oil  shale  is  contained  in  the  Parachute  Creek  member 
and  can  be  traced  over  large  portions  of  the  basin.   The 
Green  River  Formation  is  considered  to  be  middle  Eocene  in 
age  (Hail,  1974). 

The  sandy  Uinta  Formation  overlies  the  Parachute  Creek 
member  of  the  Green  River  Formation  and  is  the  dominant  unit 
in  the  tract  C-a  area.   The  Uinta  Formation  was  renamed  from 
the  Evacuation  Creek  member  (Cashion,  1974)  and  consists 
primarily  of  brown  sands  but  also  contains  lenses  of  material 
that  resembles  portions  of  the  Parachute  Creek  member. 

The  lacustrine  sediments  in  the  Piceance  Basin  were 
deposited  in  varying  stages.   As  Lake  Uinta  became  more 
alkaline,  stratification  of  the  water  developed.   As  this 
chemical  stratification  continued  the  pH  of  the  lower  level 
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Figure   2.      Plan  Map  Showing   Location   of 
Core  Holes   Used   in  the   Stratlgraphlc 
Cross  Section.    (Rohler. 197^) 
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Figure    3.      Cross   Section   of  Lithology    in    the 
Central   Portion   of   the    Piceance    Basin.    ( Sohler, 1 974 ) 
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waters  rose  until  the  lake's  lower  area  was  a  solution  of 
sodium  bicarbonate.   Organic  material  was  deposited  and 
preserved  along  with  evaporite  minerals  such  as  halite  and 
nahcolite  until  clastic  deposition  resumed  and  the  Uinta 
sediments  filled  the  basin  (Smith,  1974). 

A  basic  understanding  of  these  processes  allows  pre- 
dictions to  be  made  about  the  materials  of  environmental 
concern.   The  black  shales  in  the  Green  River  Formation  were 
deposited  in  unusually  high  pH  conditions  whereas  most 
organic  shales  are  deposited  in  acidic  conditions.   This 
explains  the  fact  that  elements  such  as  Mo,  B,  Se,  and  F 
are  elevated  in  concentration.   These  elements  generally 
become  more  mobile  in  alkaline  conditions  and  are  toxic  to 
plants  or  animals,  thus  they  are  of  prime  environmental 
concern. 

Previous  work  in  the  Piceance  Basin  (Ringrose  and 
Klusman,  1976b),  indicates  that  these  elements  are  indeed 
elevated  in  concentration  and  may  pose  pollution  hazards 
during  oil  shale  development.   Table  1  shows  some  values  of 
trace  elements  in  rocks  and  soils  from  the  earlier  studies 
and  compares  them  with  crustal  averages.   The  whole  rock 
analyses  were  determined  by  x-ray  fluorescence.   Many  of  the 
elements  vary  greatly  over  localized  areas  making  inter- 
pretations difficult  but  Zn,  Cu,  Li,  Fe,  and  Be  displayed 
significant  regional  trends  with  highest  concentrations 


Table  1 .   Trace  Elements  In  Rocks  and 
Soils  of  the  Piceance  Basin  versus 
Crustal  Averages. 


(-3 

o 
o 

Element 

Crustal 
Avg. 

Shale1 

Avg. 

Uinta2 
Formati 

on 

Parachute 
Creek 

5 

Piceance 

Basin 

Formati 

on 

Soils 

Mln.   Avg. 

Max. 

Min.   Avg. 

Max. 

Mo 

1.5 

2 

2      3 

6 

3     17 

40 

5.3 

As 

1.8 

6.6 

5     10 

21 

9     33 

110 

6.4 

Se 

0.05 

0.60 

2      2 

3 

2      2 

2 

0.28 

Zn 

70 

80 

19     39 

52 

29     47 

95 

80 

Hg 

0.08 

0.4 

0.041 

Li 

20 

60 

34 

B 

10 

100 

61 

Cd 

0.2 

0.3 

= 

1  Krauskopf,  1967 

2  Whole  Rock  Analysis  by  X-ray  FluoreBence  Klusman  and  Ringrose,  1977 

5  Ringrose  et  al,  1976b 
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observed  in  the  southern  part  of  the  basin  (Ringrose  et  al, 
1976b) . 

Both  the  Parachute  Creek  member  and  Uinta  Formation 
outcrop  in  the  vicinity  of  tract  C-a  so  a  reconnaissance 
geologic  map  (Figure  4)  of  the  study  area  was  prepared  by 
the  author  to  separate  the  sample  localities  according  to 
the  geologic  unit  from  which  they  were  derived.   The  defini- 
tion of  these  two  populations  is  difficult  because  of  the 
interfingering  of  tongues  of  Parachute  Creek  member  in  the 
Uinta  Formation.   Therefore  the  distinction  rests  upon  the 
variation  in  color  of  the  soils.   In  all  cases  the  Parachute 
Creek  soils  are  white-gray  in  color  and  in  most  instances 
the  Uinta  Formation  soils  are  brown. 
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SAMPLE  DESIGN  AND  COLLECTION 

The  sampled  area  is  eight  by  six  miles  incorporating 
tract  C-a  and  downwind  areas.   Sample  localities  were  placed 
at  half  mile  intervals  in  a  grid  pattern  within  the  area 
(Figure  5) . 

In  addition,  thirty-two  analysis  of  variance  samples 
were  selected  from  four  randomly  chosen  sections  (Figure  6) . 
Subsampling  in  each  of  the  four  sections  is  also  randomized. 
The  grid  spacing  was  determined  as  a  result  of  previous 
studies  (Ringrose,  Klusman,  1976b)  in  which  a  significant 
amount  of  the  unexplained  variance  of  trace  elements  was 
contained  at  this  scale. 

At  each  location  soil  and  three  plant  species  were 
compositly  sampled  over  an  area  of  fifty  to  one  hundred 
square  meters  depending  on  plant  availability.   Composite 
sampling  was  done  to  reduce  the  localized  (0-10m)  variance 
evident  from  previous  studies.   Plants  collected  include  big 
sage  (Artemisia  tridentata-tridentata  and  subspecies 
Artemisia  tridentata-wyomingensis) ,  Indian  rice  grass 
(Oryzopsis  hymenoides) ,  and  Western  wheat  grass  (Agropyron 
smithii) .   The  choice  of  these  particular  plants  was  based 
on  suggestions  from  range  specialists  at  Colorado  State 
University,  Area  Oil  Shale  Supervisors  Office  and  Meeker 
offices  of  the  Bureau  of  Land  Management  and  Soil  Conservation 
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Figure  5.   Plot  of  Grid  Sample  Locations 
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Service  because  of  their  availability  and  importance  as 
forage  by  cattle  and  native  herbivores,  and  potential  for 
use  in  revegetation  of  spent  shale  dumps. 

The  soil  samples  are  from  the  A  horizon  and  were  sieved 
in  the  field  with  a  4  mesh  stainless  steel  sieve  into  a 
plastic  beaker.   The  soils  were  transported  in  paper  bags 
as  was  the  plant  material. 

Stream  sediments  were  also  collected  in  most  of  the 

streams  and  gulches  up  and  downstream  from  the  grid  area  as 

well  as  within  the  area.   Stream  sediments  were  co.mposited 

over  at  least  10m  and  collected  at  one-half  mile  intervals. 
All  the  samples  were  collected  in  July  and  August, 

1976. 
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CHEMICAL  ANALYSIS 

Analytical  determinations  were  performed  by  the  Project 
Central  Analytical  Facility  at  the  University  of  Colorado 
and  also  at  Colorado  School  of  Mines.   The  elements  B,  As, 
Mo,  F,  and  Se  cannot  be  easily  analyzed  at  present  at  the 
Chemistry-Geochemistry  Department  of  Colorado  School  of 
Mines.   Project  Central  Analytical  Facility  has  completed 
analysis  only  of  B  and  Mo  in  the  soils  and  plants  because 
of  soil  matrix  interferences  that  hinder  the  analysis  for  the 
other  elements.   At  C.S.M.,  Cd,  Zn,  Hg,  and  Li  were  analyzed 
in  soils  as  elements  of  secondary  environmental  and  geo- 
chemical  interest.   The  Cd  analyses  were  abandoned  because 
of  the  strong  suspicion  of  chemical  interference  from  Ca. 
The  Ca  causes  higher  values  of  Cd  to  be  observed.   There  is 
a  strong  association  between  Zn  and  Cd  in  most  geochemical 
environments.   The  Zn  values  are  near  crustal  average  so  it 
is  expected  that  the  Cd  should  also  be  near  crustal  average. 
The  Cd  values  are  high  and  positive  correlations  exist 
between  Cd  and  Ca.   Since  Ca  is  common  in  the  marlstones  in 
the  study  area  the  Cd  results  found  were  probably  erroneous. 

All  the  elements  were  determined  by  direct  atomic 
absorption  methods  except  for  Hg  which  used  a  flameless 
atomic  absorption  method  and  Mo  in  plants  which  is  done 
colorimetrically.   The  analytical  methods  and  sample  prepara- 
tions are  described  in  detail  in  Appendix  I. 
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STATISTICAL  ANALYSIS 

Analysis  of  Variance  Design 

The  definitions  of  sample  localities  and  factors  are  as 

used  in  papers  by  Miesch  (1976a,  1976b).   The  samples  were 

collected  as  in  Figure  7.   The  sampling  was  designed  to 

analyze  the  variability  shown  in  five  separate  levels  (Figure 

8).   The  levels  are:   variance  at  the  greater  than  1.6  km 

level,  variance  at  the  0.4  km  level,  variance  at  the  0 . 1  km 

level,  variance  between  samples  at  the  50  m  level,  and  error 

in  the  chemical  analysis.   The  model  is  defined  as: 

X..,,   =  u  +  a.  +  b..  +  c.  .,  +  d.  .,  ,  +  e.  ., 
ljklm        i    ij     iik    ljkl    ijklm 

where  u  is  the  mean  of  all  the  nested  samples,  a  is  the 

regional  component  (>1.6  km),  b  is  the  0.4  km  segment,  c  is 

the  0.1  km  segment,  d  is  the  50  m  segment,  and  e  is  the 

analysis.   The  population  variance  is  divided  as  follows: 

2     2     2     2     2     2 
x    a    b    c    d    e 

and  is  calculated  as  the  sample  variance: 

S2  =  S2  +  S2  +  S2  +  S2  +  s2 
x    a    b    c    d    e 

The  complete  analysis  of  variance  data  are  listed  in 

Tables  11-14  in  Appendix  II.   In  Table  2  the  geometric  means 

and  deviations  and  variance  ratios  are  listed.   The  methods 

used  to  calculate  the  geometric  means  and  deviations  are 

given  in  Appendix  II.   In  the  interpretation  of  the  analysis 
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of  variance  design,  ratios  of  the  variances  at  different 
levels  of  the  sampling  model  are  examined.   The  variance 
ratio  (V)  is  the  ratio  of  the  variance  between  sections 
(>1.6  km)  and  the  sum  of  the  variance  at  the  other  levels 
of  the  model  (0.4  km,  0.1  km,  50  m,  and  analytical  error). 
The  larger  the  variance  ratio  the  greater  the  probability 
that  the  component  in  question  displays  a  regional  compon- 
ent of  variation.   The  variance  mean  ratio  (V  )  is  similar 
to  the  variance  ratio  but  incorporates  the  nested  design  of 
the  sampling  model  used  (see  Appendix  II).   The  maximum 
acceptable  error  variance  for  a  balanced  sampling  design 
(Er)  and  the  error  variance  for  a  hierarchial  design  (E  ) 
are  used  to  determine  the  stability  of  a  geochemical  map 
(Miesch,  1976a) .   If  the  variance  observed  in  a  nested 
design  (Es)  is  less  than  the  maximum  permissible  error  var- 
iance for  a  balanced  design  (E  ) ,  then  the  model  produces 
an  accurate  representation  of  the  variance  shown  by  the 
particular  component  (Miesch,  1976b).   This  also  supports 
the  stabilities  of  the  estimates  of  the  means  and  suggests 
there  is  a  significant  regional  component  of  variance 
(>1.6  km)  for  As,  Li,  Mo,  and  organic  carbon  in  soil  and  B 
and  Mo  in  sage.   If  the  variance  ratio  is  zero  it  implies 
there  is  no  regional  component  of  variance  and  the  maps  for 
these  components  are  not  stable.   This  is  the  case  for  Hg, 
B,  and  pH  in  soil  and  Hg,  Zn,  B,  and  Mo  in  wheat  grass  and 
also  Zn  in  rice  grass  (Table  2) . 
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Table  2.   Geometric  Means  and  Deviations  and 
Variance  Ratio  Results  for  Tract  C-a  Data. 


ELEMENT    95^ 

AND  RANGE 
MEDIA 

WHEATGRASS 

Hg(ppb)         7.5-70 

Zn(pmn)        5.1-18 

B(ppm)  7.7-28 

Mo(ppm)        0.54-2.** 


GEOMETRIC      GEOMETRIC      VARIANCE**  VARIANCE      N, 


MEAN 


DEVIATION 


23 

1.75 

9.6 

1.38 

15 

1.38 

1.2 

1.45 

RATIO 


MEAN 
RATIO 


E. 


# 
Of 

SAMPLES 


25 
25 
25 
25 


H9 
IV) 

o 

© 


RICEGRASS 
Ha;(ppb)         11-52 
Zn(ppm)      0.49-18 
B(ppm)  3-32 

Mo(ppm)      0.48-2.3 


24 

1.47 

2.6 

2.46 

10 

1.81 

1.1 

1.48 

0.05     2.97   6o  0.00059  0.00061  32 

32 

0.02       2.96    45   0.0004    0.0001    32 
0.24      3.95    4  0.0068   0.0017   32 


Arithmetic  mean  and  deviation  for  pH  (r»H  Is  a  loo;  measurement) 

** 

If  the  estimated  variance  between  sections  Is  zero  V  ,  N  ,  E  ,  E  ore  not 

m   r   r   8 

calculated . 


(VI 
(V) 


Table  2.   Continued  . 


ELEMENT 

AND 

MEDIA 

SOILS 

Hg(ppb) 

Zn(ppra) 

Li(ppm) 

B(ppm) 


95* 
RANGE 


30-58 
39-123 
9-40 

78-195 


GEOMETRIC      GEOMETRIC      VARIANCE **VARIANCE      N 


Mo(ppra)      0.46-4.7 

As(ppra)        4-18 

Organic 

Carbon  %     0.46-2.9 


pH 

sIge 

Zn(ppm) 

B(ppm) 

Mo(ppm) 


7.1-8.5 


MEAN 

42 

70 

20 

123 

1.6 

8.7 

1.2 

7.8* 


DEVIATION 


RATIO 


0.43-11.         2.2 

21.5-^2.6     30.3 
C. 25-1. 7        O.65 


1.18 

1.33 
1.43 

1.26 

1.70 

1.43 

1.58 
0.353* 

2.26 
1.19 
1.63 


0.02 
0.41 

0.23 
1.36 

1.03 


0.02 
0.46 
0.?2 


MEAN 
RATIO 


4.17 
4.32 

4.28 
4.40 

3.96 


2.96 
3.74 
3.90 


45 
3 

4 
2 


E 


OF 
SAMPLES 


253 
0.019        0.0002      253 

0.011       0.003       253 

253 

0.011      0.003      253 

0.0071   0.0048   32 
0.013   0.010   2^3 

253 


45  0.004  0.001  32 
2  0.0024  0.0006  249 
5    0.0089   0.0026   243 


rv> 
O 
O 
v-D 
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The  lack  of  stability  for  Hg  and  Zn  in  the  soils  and 
plants  is  probably  due  to  a  high  component  of  variance 
caused  by  error  in  the  methods  of  chemical  analysis.   Many 
of  the  elements  analyzed  in  wheat  grass  show  little  stabil- 
ity because  the  power  of  the  model  is  reduced.   At  many  of 
the  analysis  of  variance  sample  sites  (also  many  of  the 
grid  localities)  wheat  grass  was  not  available  for  collec- 
tion, thus  the  model  is  less  likely  to  produce  non-zero 
estimates  of  variance  at  the  highest  level  of  the  model. 

All  parameters  that  have  non-zero  regional  variance 
components  show  enough  variation  for  geochemical  maps  to  be 
drawn.   In  most  cases  those  elements  with  zero  regional 
components  have  a  high  analytical  error  component.   There- 
fore the  choice  of  sampling  interval  and  grid  spacing  was 
sufficient  to  describe  most  of  the  variance  in  this  area. 

An  important  step  in  statistical  analysis  is  to  examine 
the  distributions.   The  normal  convention  with  trace  element 
data  is  to  transform  the  parts  per  million  by  a  logarithmic 
(base  10)  conversion.   Figure  9  displays  histograms  of  Mo  in 
the  soil  samples,  first  in  ppm,  then  in  log,Q  ppm.   A 
logarithmic  transformation  noticeably  improves  the  distribu- 
tion to  a  more  nearly  normal  distribution.   The  other  ele- 
ments studied  display  similar  improvements  in  their  distribu- 
tions. 
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The  next  step  is  to  compare  the  means  and  deviations  of 
the  data.   Data  collected  previously  in  the  Green  River  Basin 
of  Wyoming  and  Piceance  Basin  is  shown  in  Table  3. 

The  data  in  this  survey  fall  within  expected  ranges  for 
all  of  these  elements  except  possibly  the  B  in  sage  and  Zn 
in  sage.   The  low  mean  of  Zn  (2.2  ppm)  in  sage  is  probably 
due  to  analytical  errors  caused  by  poor  digestions  of  sage 
with  perchloric  acid  (see  Appendix  I).   The  B  results  (30 
ppm)  are  higher  than  expected  and  will  be  discussed  later. 
No  data  could  be  found  with  rice  grass  or  wheat  grass  from 
similar  geographic  areas  for  comparison  but  the  grasses  are 
elevated  in  Mo  with  respect  to  the  sage  at  the  same  sample 
site.   The  toxic  level  in  forage  for  cattle  is  5-6  ppm  and 
the  grasses  average  from  1.1  to  1.2  ppm,  so  no  immediate 
problem  exists  (Gough,  1976). 

Table  4  defines  the  correlation  matrix  for  the  compon- 
ents analyzed  and  the  correlations  are  as  expected.   Organic 
carbon  and  pH  are  negatively  correlated  and  Mo  in  soil  and 
Mo  in  sage  are  strongly  correlated.   Li,  Mo,  and  B  also 
exhibit  a  strong  positive  relationship  in  the  soil  samples. 
The  only  inconsistency  is  the  weakly  positive  relationship 
between  boron  in  soils  and  boron  in  sage.   This  is  discussed 
in  the  hypothesis  testing  section. 

The  samples  came  from  more  than  one  lithology  and  there- 
fore more  than  one  population.   This  causes  problems  in  the 
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Table  3.   Trace  Elements  in  Soils 
of  the  Piceance  Fasin  and  Sagebrush 
of  the  Green  River  Basin,  (ppra) 
(TT.S,  Geological  Survey  Open  File 
Report  76-729,  Appendix  Til). 


Element 

S 

DilS 

Hg 

0 

.041 

Zn 

80 

As 

19 

B 

61 

Mo 

5.3 

Li 

34 

Cd 

„ 

Sagebrush(dry  weight) 

0.034 

28 

0.64 

13 

0.70 

1.4 

0.34 


Table  4.  Correlation  Matrix 
of  Tract  C-a  Vicinity  Soil 
Samples. 


1-3 

o 

o 

U3 


Variable 

Hg 

1.000 

Zn 

0.026 

1.000 

Li 

0.092 

0.137 

1.000 

Org-C 

0.140 

0.055 

0.395 

1.000 

PH 

-0.287 

0.101 

0.077 

-0.143 

1.000 

Sage  MO 

0.130 

0.011 

0.378 

-0.004 

0.075 

1.000 

Sage   B 

0.083 

0.002 

-0.052 

-0.114 

-0.011 

0.037 

1.000 

B 

0.007 

0.048 

0.274 

0.171 

0.089 

0.161 

0.053 

Mo 

0.119 

0.100 

0.605 

0.143 

0.008 

0.419 

-0.061 

Hg 

Zn 

Li 

Org-C 

PH 

Mo 
Sage 

B 
Sage 

1.000 

0.2156         1.000 

B  Mo 
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interpretation  of  the  analysis  of  variance  data  because 
the  techniques  are  developed  to  consider  only  single  popu- 
lations.  It  was  suspected  that  the  Parachute  Creek  member 
samples  are  different  from  the  Uinta  Formation  samples. 
Differences  in  trace  element  concentrations  are  evident 
from  whole  rock  analyses  (Table  1).   This  would  presumably 
affect  the  regional  component  of  variance  in  soils  because 
the  Parachute  Creek  member  outcrops  only  in  the  western 
segment  of  the  study  area.   The  reconnaissance  geological 
map  (Figure  4)  allowed  the  separation  of  samples  according 
to  lithology.   Hypothesis  testing  can  be  employed  to  test 
if  the  differences  in  means  of  the  samples  on  the  different 
formations  are  significant.   Checks  on  the  variances  of  both 
populations  revealed  no  significant  differences  so  the 
Student's  t  test  can  be  used. 

Hypothesis  Testing 

The  results  of  t  tests  are  given  in  Table  5.   Table  5a 
shows  that  all  the  parameters  except  Hg  and  Zn  in  soils  and 
the  pH  of  the  soils  show  significant  differences  between 
the  Parachute  Creek  member  and  Uinta  Formation. 

Since  there  is  a  difference  between  these  two  popula- 
tions the  next  logical  question  is  if  there  are  any  trends 
in  either  population  over  the  study  area.   The  Parachute 
Creek  member  is  so  locally  distributed  that  it  is  probably 
safe  to  assume  that  there  is  little  change  in  the  concentra- 


T  2009  30 


Table  5a.   Uinta  Formation  Samples  vs  Parachute  Creek  Member 
Samples 


Element 

X  Uinta 

s 

X  Par.  Crk. 

3 

t  Value 

DF 

Hg 

45.9 

1  .40 

50.1 

1  .50 

1.12 

39 

Zn 

69.2 

1.30 

73.0 

1.35 

1.27 

42 

Li 

17.8 

1.43 

30.2 

1.48 

6.62** 

42 

B 

121.2 

1.27 

133.4 

1  .40 

2.43** 

42 

Mo 

1.44 

1.51 

2.72 

1.59 

6.14** 

42 

Mo  (sage) 

0.61 

0.60 

0.88 

0.71 

3.77** 

40 

B  (sage) 

30.0 

1.53 

31.9 

1.48 

2.46** 

41 

Org  C  (%) 

1.08 

1.58 

1.64 

1.56 

6.22** 

42 

pH 

7.85 

0.35 

7.79 

0.32 

1.00 

42 

Table  5b.  Test  Results  of  Uinta  Formation  Samples,  Group  1 

vs  Group  2. 

Element  "~  Gronpi    s      X  Grouo  2  3  t  Value  DF 

Hg  45-°     1.33       43.0  1.38  0.64  56 

TZ?  70-1     1.42       66.3  1.45  1.11  79 

i1  17-3    1.50       17.2  1.46  0.14  79 

B  122-°     1.30      121.0  1.32  0.24  79 

™>     „  1-47    1.47        1-39  1.52  0.97  79 

B  (sage)  29.2     1.45       30.4  1.43  1.56*  77 

Mo  (sage)  0.63   o.58        0.59  0.63  0.77  77 

Table  5c.  t  Test  Results,  Uinta  Formation  Samples,  Group  1 

vs  Group  5. 

Element  f  CLUQUP  1    3     I  Crmip  3  s  t  Value  DF 

HS  45.0      1#53     45.9  1.30  0.27  47 

Zn  70.1       1#42      66.4  1.40  1.25  49 

J;3-  17-3      1.50     19.6  1.43  2.89**  49 

B  (soil)  122.0      1#30    125.O  1.35  0.63  49 

Mo  (soil)  1.47     1#47      1.58  1.42  0.79  49 

B  (sage)  29.2      } # 45     31.7  1.41  2.78**  49 

Mo  (sage)  0.63     0.58      0.62  0.61  0.08  49 

Org  C  (I)  1.04     U5Q               1.22  1.47  2.13**  49 

PH  7-9      0.32       7.6  0.28  3.84**  49 


**signif icant  at  a  =  .05 
*significant  at  a  =  .10 


Groups  1-4  are  defined  in  Appendix  III  pp. 95-98. 
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Table  5d.   t  Test  Results,  Uinta  Formation  Samples,  Group  2 
vs  Group  5. 


Element 

X  Grout)  2 

s   ■  X 

Group  3 

s 

t  Value 

DF 

Hg 

42.9 

1.38 

45.9 

1.30 

0.85 

47 

Zn 

66.2 

1.45 

66.4 

1  .40 

0.03 

49 

Li 

17.2 

1.46 

19.6 

1.43 

2.62** 

49 

B 

121.0 

1.32 

125.0 

1.35 

0.68 

49 

Mo 

1.38 

1.52 

1.58 

1.42 

1.36* 

49 

B  (sage) 

30.4 

1.43 

31.7 

1.41 

1.52* 

49 

Mo  (sage) 

0.59 

0.63 

0.62 

0.61 

0.58 

49 

Org  C  (1) 

1.02 

1.53 

1.22 

1.47 

2.33** 

49 

pH 

7.9 

0.33 

7.6 

0.28 

4.87** 

49 

Table  5e.   t  Test  Results,  Uinta  Formation  Samples,  Group  3 
vs  Parachute  Creek  Member  Samples,  Group  5 

Mean  Par.        Mean,  Group 
Element   Crk.  Samples   s    3  (Unita)      3   .t  Value    DF 

Hg  50.1  1.50  45.9  1.30  0.93  39 

Zn  73.0  1.35  66.4  1.40  2.46**  42 

Li  30.2  1,48  19.6  1.43  4.49**  42 

B  133.4  1.40  125.0  1.35  1.31*  42 

Mo  2.7  1.59  1.6  1.42  4.11**  42 

B  (sage)  31.9  1.48  31.7  1.41  0.21  41 

Mo  (sage)  0.88  0.71  0.62  0.61  3.09**  40 

Org  C  {%)  1.64  1.56  1.23  1.47  3.59**  42 

pH  7.8  0.32  7.6  0.28  2.18**  42 

Table  5£.   t  Test  Results,  Sagebrush  Subspecies  Tridentata 
vs  Subspecies  Wyomingensis . 

Element    X  Tridentata   s    X  '.Vyomiriiger.sis   s     t Value   DF 

Mo  (sage)  0.64  0.48  0.67  0.54  0.54  43 

B  (sage)  30.4  1 . 50  30.3  1.47  0.11  44 

Mo  (soil)  1.52  1.53  2.07  1.54  2.73**  45 

B  (soil)  121.9  -1,25  129.2  1.38  1.38*  45 

Table  5g.   t  Test  of  Sagebrush  Subspecies  Wyomingensis  on 

Parachute  Creek  Member  vs  Subspecies  Wyomingensis 
on  Uinta  Formation. 

Wyom.  Mean  Wyom.  Mean 

Element   Parachute  Crk     3     Uinta  Fm      3       Value  DF 

Mo  (sage)  0.94  0.50  0.57  0.61  3.34**  14 

B  (sage)  30.9  1.36  29.9  1.54  0.74  15 

Mo  (soil)  3.77  1.40  1.46  1.36  4.58**  16 

B  (soil)  143.3  1.36  121.7  1.44  2.73**  16 
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tions  over  such  a  small  geographic  distance.   The  Uinta 
Formation  is  spread  throughout  the  whole  area  so  artificial 
boundaries  were  created  that  separated  the  Uinta  Formation 
samples  into  three  groups  shown  in  Figures  18  to  20 
(Appendix  III),  where  group  1  is  in  the  easternmost  part 
of  the  area  and  groups  2  and  3  proceed  westerly.   Table  5b 
shows  a  test  of  group  1  vs  group  2.   Only  boron  in  sage 
shows  a  difference  in  means  at  the  or-  .10  level.   Table  5c 
shows  the  test  between  group  1  and  group  3  on  Uinta  Formation 
(max.  geographic  separation).   In  this  case,  Li,  B  in  sage, 
organic  carbon,  and  pH  in  soils  show  significant  differences. 
In  Table  5d  group  2  vs  group  3  Uinta  Formation  samples  are 
tested.   Li  and  Mo  in  soils,  B  in  sage,  organic  carbon,  and 
pH  have  significantly  higher  means.   The  means  of  many  com- 
ponents in  the  Uinta  Formation  are  increasing  westward. 
These  trends  in  the  Uinta  are  most  strongly  controlled  by 
organic  carbon  and  pH.   Therefore  the  distinctions  are 
probably  caused  chiefly  by  current  climatic  differences, 
since  the  western  area  has  higher  elevations,  more  moisture, 
and  a  thicker  vegetation  cover. 

These  climate  variations  are  significant  enough  to 
cause  differences  within  the  Uinta  Formation  and  they  are 
further  emphasized  by  observed  differences  between  the 
Parachute  Creek  member  and  the  Uinta  Formation.   Table  5e 
shows  a  test  of  Parachute  Creek  samples  vs  Uinta  samples 
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both  in  the  same  geographic  area  (Group  3) .   All  components 
except  Hg  and  B  in  sage  show  significantly  higher  means  in 
the  Parachute  Creek  member  samples.   Therefore  the  two 
lithologies  do  have  differences  in  trace  elements  in  their 
respective  soil  profiles. 

In  order  to  test  the  appropriateness  of  the  artificial 
distinctions  of  groups  1-3  each  of  the  three  groups  were 
randomized  by  random  deletions  and  the  t  tests  were  rerun. 
This  was  done  to  minimize  the  possibility  of  single  sample 
localities  influencing  the  population  means.   In  all  of  the 
randomized  tests  no  differences  in  t  test  data  were  found. 

Another  test  was  devised  to  ascertain  whether  species 
differences  are  apparent  in  sagebrush.   The  two  subspecies 
of  sage  (wyomingensis   and  tridentata)  were  tested  and 
Table  5f  shows  that  no  significant  difference  exists  in  the 
subspecies  even  though  the  B  and  Mo  are  significantly 
different  in  the  soils  on  which  they  grew.   This  suggests 
that  these  two  subspecies  of  sage  vary  little  over  the  study 
area.   Table  5g  shows  the  means  of  the  wyomingensis  with 
respect  to  differences  in  lithology.   This  shows  that  the  B 
and  Mo  in  the  soils  of  the  Parachute  Creek  and  Uinta 
Formation  are  significantly  different  in  this  test.   However, 
only  the  Mo  shows  a  significantly  different  mean  in  the  sage. 
This,  in  part,  may  be  due  to  the  fact  that  it  is  only  a  small 
population  (df  =  15) . 
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Many  inconsistencies  exist  with  respect  to  the  B  values 
in  sage.   Boron  in  sage  has  a  high  regional  component  of 
variance  while  boron  in  soil  has  no  regional  component  of 
variance.   Boron  in  sage  has  a  negative  correlation  with 
organic  carbon  in  the  soil  but  B  in  soil  has  a  positive 
correlation  with  organic  carbon  in  soil.   This  situation 
spawns  speculation  about  B  in  sage.   These  inconsistencies 
could  be  due  to  a  biochemical  influence.   Table  3  shows 
that  other  data  from  the  similar  areas  list  B  in  the  soil 
as  61  ppm  (GM)  and  in  sage  as  13  ppm  (highest  expected 
value).   Since  in  the  vicinity  of  tract  C-a,  both  the  soils 
and  sage  are  much  higher  it  might  be  possible  that  the  sage 
is  exhibiting  some  sort  of  leveling  effect.   As  the  B  in 
soil  increases  the  B  in  sage  increases,  but  at  a  slower  rate, 
This  could  be  why  B  in  sage  exhibits  a  regional  trend 
whereas  B  in  soil  has  none. 

The  slower  increase  in  boron  in  sage  may  be  dampening 
the  local  variance  exhibited  by  B  in  soil.   This  could  also 
be  the  reason  why  the  B  in  soil  and  sage  do  not  correlate 
well.   It  must  be  emphasized  that  this  is  speculation  and 
additional  studies  would  be  necessary  to  show  any  of  these 
relationships  conclusively. 
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Table  6  .   Results  of  Multiple  Regression  with 
B  as  the  Dependent  Variable. 


Sample  Size 

Dependent  Variable 

Independent  Variables 

Coefficient  of 
Determination 

Multiple  Corr 
Coefficient 

Estimated  Constant 
Term 

Standard  Error  of 
Estimate 

Analysis  of  Variance 
for  the  Regression 
Source  of  Variation 

Regression 

Residuals 

Total 

REGRESSION 
VAR.     COEFFICIENT 

Hg  0.2364452E-01 

Li  0.1186452 

Org  C  0.364379BE-02 

pH  0.3552470E-01 

B  -0.4397217E-02 

MoSo  0.2825871E-02 


253 

B 

Hg,  Zn, 

0.09282 

0.30466 

1.7479070 

0.10072519 


Li,  Org,  pH,  Mo 


DF 

6 

246 
252 


S.  SQ. 

0.255365 

2.49581 

2.75117 


M.S.        F      PROB 

0.425608E-01   4.195   0.0005 
0.101456E-01 


S.E.  OF 
REG.  COEF, 


F- VALUE 
OF  (1.246) 


0 
0 
0 
0 
0 
0 


3833E-01 
3833E-01 
4391E-01 
2311E-01 
7901E-01 
4241E-01 


0.3805 


2. 

0 

2, 

0, 

0 


02 


863 

6888E- 

363 

3831E-02 

4440E-02 


PROB 

5379 
0919 
9339 
1256 
9507 
9469 


CORR.  COEFF. 
WITH  B 


0 
0 
0 
0 
0 
0 


0262 
1637 
0553 
1014 
0480 
1003 


T  2009 


36 


Table  7  .   Results  of  Multiple  Regression  with 
Zn  as  the  Dependent  Variable. 


Sample  Size 

Dependent  Variable 

Independent  Variables 

Coefficient  of 
Determination 

Multiple  Corr 
Coeff 

Estimated  Constant 
Term 

Standard  Error  of 
Estimate 

Analysis  of  Variance 
for  the  Regression 
Source  of  Variation 

Regression 

Residuals 

Total 


VAR. 

Hg 

Zn 

Li 

Org  ( 

PH 

M0S0 


REGRESSION 
COEFFICIENT 

0.4432343E-02 

0.3184751E-02 

0.1105455 

0.5443270E-01 

0.2553415E-01 

0.4465163E-01 


253 

Zn 

Hg,  Li,  Org,  pH,  B,  Mo 

0.03634 

2.19062 

1.3832824 

0.12481431 


DF 

6 
246 
252 

S.E, 
REG. 


S.  SQ. 

0.144502 

3.83234 

3.97684 


M.S.  F      PROB 

0.240837E-01   1.546   0.1638 
0.155786E-01 


OF      F-VALUE 

COEF.  OF  (1.246)  PROB 

0.3096E-01  0.2050E-01  0.8863 

0.5145E-01  0.3831E-02  0.9507 

0.5648E-01  3.832  0.0514 

0.3526E-01  2.383  0.1240 

0.1867E-01  1.871  0.1727 

0.3410E-01  1.714  0.1917 


CORR.  COEFF, 
WITH  Zn 


0 
0 
0 
0 
0 
0 


0072 
0480 
2736 
1710 
0896 
2156 
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Table  8  .   Results  of  Multiple  Regression  with 

Mo  as  the  Dependent  Variable. 


Sample  Size 

Dependent  Variable 

Independent  Variables 

Coefficient  of 
Determination 

Multiple  Corr 
Coeff . 

Estimated  Constant 
Term 

Standard  Error  of 
Estimate 

Analysis  of  Variance 
for  the  Regression 
Source  of  Variation 

Regression 

Residuals 

Total 


253 

Mo 

Hg,  Zn,  Li,  Org,  pH,  B 

0.38820 

0.62305 

-1.2369728 

0.18765236 


DF 

6 

246 


S.  SQ. 

5.49647 
8.66250 


M.S. 

0.916078 
0.352134E- 


F 
26.02 


PROB 
0.0000 


01 


REGRESSION 


252  14.1590 
S.E.  OF 


F- VALUE 


VAR 

Hg 

Zn 

Li 

Org 

pH 

B 


COEFFICIENT    REG  COEF.    OF  (1.246) 


0.7158407E-01  0 

0.6387511E-02  0 

0.9662530  0 

0.1604181  0 

0.3127205E-01  0 

0.1549777  0 


5749E-01 
9586E-01 
8630E-01 
6521E-01 
3486E-01 
1184 


1.550  0 

0.4440E-02  0 

125.4  0 

6.051  0 

0.8049  0 

1.714  0 


PROB 

2143 
9469 
0000 
0146 
3705 
1917 


CORR.  COEFF, 
WITH  Mo 


0 
0 
0 
0 
0 
0 


1192 
1003 
6047 
1432 
0051 
2156 
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Table  9  .   Results  of  Multiple  Regression  with 
Li  as  the  Dependent  Variable. 


Sample  Size 

Dependent  Variable 

Independent  Variables 

Coefficient  of 
Determination 

Multiple  Corr. 
Coeff . 

Estimated  Constant 
Term 

Standard  Error  of 
Estimate 

Analysis  of  Variance 
for  the  Regression 
Source  of  Variation 

Regression 

Residuals 

Total 


VAR. 

Hg 
Zn 
Org 

PH 
B 

Mo 


REGRESSION 
COEFFICIENT 


253 

Li 

Hg,  Zn,  Org,  pH,  B, 

0.49132 

0.70094 

0.36804895 

0.11283737 


Mo 


DF 

6 
246 
252 

S.E. 
REG. 


S.  SQ. 

3.02526 
3.13214 
6.15740 


M.S. 

0.504210 
0.127323E-01 


F 
39.60 


PROB 
0.0000 


OF      F-VALUE 
COEF.   OF  (1.246) 


0.9696773E-01  0 

0.2418713  0 

0.4486958E-01  0 

0.1387339  0 

0.3493725  0 


5731E-01 
3657E-01 
2080E-01 
7087E-01 
3120E-01 


2.863 
43.73 
4.654 
3.832 
105.4 


PROB 


0.8781140E-02   0.3467E-01   0.6413E-01   0.8003 


0 
0 
0 
0 
0 


0919 
0000 
0319 
0514 
0000 


CORR  COEF 
WITH  Li 


0 
0 
0 
0 
0 
0 


0919 
1637 
3955 
0766 
2736 
6047 
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Figure  12.   Trend  Map  of  Organic  Carbon  in 
all  Soil  Samples.  (degree=3) 
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Table  10.   Trend  Surface  Analysis 
Statistics. 


Element 


Trend 
Degree 


F-value 


LI 


10.9*(10,2^2)     0.^6 


LI (Uinta) 


7.3*(10.199)     0.52 


Ko 


6.2*(10,2^2)     0.-35 


Mo(Uinta) 


2.8*(10,199)     0.36 


B(sage) 


if. 17(2, 237)      0.20 


Organic 
Carbon 


6.1*(10,2^2)      0.^5 


Mo (sage) 


2.o*(6;236)     0.22 


*  Significance  at  <*  =0.05 

Multiple  Correlation  Coefficient  is  equal  to  the  square 
root  of  the  cercent  sum  of  sauares  explained  by  the  re- 
aression. 


c 
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Other  Data  Reduction  Techniques 

Multiple  regressions  (Tables  6-9)  were  also  run  with 
the  important  components  as  dependent  and  alternately  as 
independent  variables.   These  serve  to  further  illustrate 
and  clarify  the  correlations  between  the  variables.   The 
only  statistically  significant  regressions  were  for  the 
dependent  variables  Li,  B,  and  Mo  in  the  soils  with  the 
greatest  influence  being  organic  carbon,  especially  in  the 
Li  regression.   The  independent  variables  from  the  soil  data 
were  B,  Hg,  Zn,  Li,  organic  carbon,  and  pH. 

Trend  surface  analyses  pictorially  display  (Figures  10- 
17)  the  analysis  of  variance  results  on  the  sample  grid. 
For  elements  with  no  regional  variance  the  trend  surfaces 
are  random  noise.   Contour  intervals  were  chosen  by  dividing 
the  range  of  values  for  each  component  by  an  integer,  whose 
value  depends  on  the  regional  variation  shown.   The  greater 
the  range  and  regional  component,  the  more  contour  intervals 
can  be  shown.   The  trend  maps  for  Li,  Mo,  and  organic  carbon 
in  soils  show  statistically  significant  trends  when  the 
lithology  differences  are  not  taken  into  account.   When  only 
the  samples  taken  on  the  Uinta  Formation  are  mapped,  Li,  Mo, 
and  organic  carbon  in  the  soils  still  show  significant  trends 
Thus  a  regional  trend  is  displayed  in  the  Uinta  Formation 
from  east  to  west  with  increasing  concentrations  toward  the 
west.   This  trend  is  enhanced  by  the  lithologic  control  of 
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the  Parachute  Creek  member.   The  Mo  and  B  in  the  sage  do  not 
exhibit  95%  confidence  in  their  trend  surfaces  but  they  are 
very  similar  to  the  soil  trends  (Table  10). 

Previous  work  in  the  Piceance  Basin  (Ringrose  et  al, 
1976b)  suggested  that  the  regional  trends  of  Li  and  Zn  were 
north  to  south  with  higher  concentrations  in  the  south  and 
west.  This  study  indicates  that  in  the  tract  C-a  area  the 
trend  has  a  significant  east-west  component.  The  trend  of 
Zn  was  statistically  insignificant,  thus  it  seems  that  for 
Zn  the  trend  is  broader  and  encompasses  a  larger  area  than 
was  sampled  in  this  study. 
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CONCLUSIONS 

A  geochemical  baseline  study  of  soils  and  plants  from 
oil  shale  tract  C-a  and  vicinity  shows  that  the  elements 
As  and  B  are  elevated  with  respect  to  normal  soils.   These 
elements  may  pose  environmental  hazards  as  the  development 
of  the  oil  shale  industry  proceeds.   The  elements  Zn,  Li, 
Mo,  and  Hg  are  present  in  average  concentrations.   The 
•geometric  means  of  Li  and  Mo  show  no  significant  difference 
from  crustal  averages  but  when  only  the  samples  from  the 
Parachute  Creek  member  are  considered,  Li  and  Mo  are  ele- 
vated.  These  elements  are  potentially  toxic  during  the 
storage  of  raw  shale  and  as  airborne  pollutants  from  retort- 
ing processes.   The  analysis  of  plant  materials  (big  sage, 
Western  wheatgrass,  Indian  ricegrass)  show  correlations  with 
the  soils  from  which  they  are  derived.   Mo  is  not  elevated 
in  sage  compared  to  other  areas  but  the  B  is  substantially 
higher  in  soils  and  sage  than  the  Piceance  Basin  as  a  whole. 

Statistical  analysis  of  means  show  that  there  is  a 
difference  between  the  soils  and  plants  on  the  Uinta  Forma- 
tion as  opposed  to  the  Parachute  Creek  member.   Most  of  the 
surficial  materials  show  significantly  higher  values  in  the 
Parachute  Creek  member.   Hypothesis  tests  also  indicate  that 
there  is  no  significant  difference  in  the  sage  subspecies 
Artemisia  tridentata  wyomingensis  and  tridentata. 
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The  limited  data  from  the  rice  grass  and  wheat  grass 
show  enrichment  in  Mo  and  B.   Future  studies  would  benefit 
by  additional  analyses  of  these  important  forage  materials. 

Regional  trends  with  increasing  concentrations  from 
east  to  west  are  displayed  within  the  study  area  for  Li,  Mo, 
and  organic  carbon  in  soils  and  B  and  Mo  in  sage.   These 
trends  are  evident  in  the  soils  developed  on  the  Uinta 
Formation.   The  regional  trend  is  enhanced  by  lithologic 
variations.   Variations  in  the  lacustrine  environment  of 
ancient  Lake  Uinta  are  likely  the  cause  for  these  regional 
variations,  since  the  study  area  is  not  located  in  the 
center  of  the  Piceance  Basin. 

The  sampling  design  was  sufficient  to  describe  most  of 
the  important  variation  in  the  environment.   Composite 
sampling  decreased  the  low  level  variance,  but  for  some 
components  the  local  variance  is  still  dominant  suggesting 
that  over  short  geographic  distances  the  area  may  seem 
heterogeneous  but  in  most  cases  the  actual  deviations  are 
minimal. 
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APPENDIX  I 

Sample  Preparation 

Each  soil  sample  consisted  of  approximately  200  grams. 
The  whole  sample  was  repeatedly  split  with  a  riffle  splitter 
until  about  10  grams  remained.   The  splitting  was  done  to 
avoid  analytical  errors  caused  by  stratification  of  soil 
components  during  transport.   The  10  gram  split  was  then 
ground  for  about  8  minutes  in  a  tungsten  carbide  vial  with 
a  leucite  sleeve  by  a  Spex  mixer  mill. 

Plant  samples  were  allowed  to  dry  at  ambient  laboratory 
temperatures.   The  wheat  grass  and  rice  grass  were  washed 
with  distilled  water  in  an  ultrasonic  bath.   The  grass 
samples  again  were  allowed  to  dry  for  three  days.   The 
washing  procedure  was  carried  out  to  remove  soil  from  the 
grass  caused  by  rain  splash.   The  sage  samples  were  not 
washed  because  only  new  growth  was  sampled  from  upper  parts 
of  the  bush  where  rain  splash  and  wind  blown  soil  do  not 
appear  to  be  significant.   All  the  plant  samples  were  ground 
in  a  commercial  blender.   The  sage  was  ground  finely  while 
the  grasses  remained  somewhat  coarser.   After  grinding  the 
samples  were  dried  at  70°C  for  24  hours. 

Analyses  were  performed  on  all  sage  samples  but  only  on 
analysis  of  variance  wheat  grass  and  rice  grass  because 
collection  ratios  of  the  grasses  tended  to  be  low,  especially 
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for  the  wheat  grass.   All  samples  were  randomly  chosen  for 
analysis  to  avoid  any  non-random  analytical  errors.   Random- 
ization is  necessary  because  non-random  analytical  errors 
could  contribute  to  regional  variations  while  analyzing 
samples  from  a  specific  geographic  area. 

To  analyze  for  Hg,  Zn,  and  Li  in  the  soil  and  plant 
samples  approximately  1  g  of  sample  was  digested  with  5:1 
perchloric-nitric  acid  for  3  hours  at  a  constant  temperature 
of  95  C.   Each  sample  was  diluted  to  50  mis  and  mixed 
thoroughly.   The  Hg  analysis  was  carried  out  immediately 
after  cooling  to  avoid  losses  of  this  volatile  component. 
Zn  and  Li  were  analyzed  from  the  same  digestion  on  the 
following  day.   The  flameless  Hg  method  (Hatch  and  Ott, 
1968)  has  a  sensitivity  of  about  20  ppb  in  1  g  of  soil.   A 
dilute  solution  of  stannous  chloride  is  added  to  an  aliquot 
of  the  sample  to  reduce  the  mercury  to  the  elemental  vapor 
state.   The  vapor  is  then  forced  through  a  quartz  tube 
positioned  in  the  beam  path.   The  absorption  is  recorded  on 
a  strip  chart  and  compared  to  standards.   All  atomic 
absorption  analyses  were  performed  on  a  Perkin-Elmer  Model 
303  atomic  absorption  spectrometer. 

The  Zn  and  Li  analyses  are  performed  by  normal  atomic 
absorption  methods.   It  should  be  emphasized  that  all 
determinations  are  total  concentrations  from  bulk,  not 
laboratory  sieved  fractions. 
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pH  and  Organic  Carbon  Analysis 

The  pH  was  determined  by  preparing  a  saturation  paste 
using  the  bulk  soil  sample  and  distilled  water.   A  combina- 
tion pH  electrode  was  immersed  in  the  sample  and  a  reading 
is  made.   Problems  with  this  method  arise  because  soils 
with  varying  grain  size  and  clay  content  may  require  differ- 
ent amounts  of  water  necessary  to  form  a  paste.   This 
causes  fluctuations  in  the  observed  pH. 

Organic  carbon  was  determined  by  the  Walkley-Black 
Method  (American  Society  of  Agronomy,  1965).   Ground  and 
accurately  weighed  samples  of  soil  were  suspended  in  10  ml 
of  1.000N  K2Cr07.   Twenty  ml  of  concentrated  H2S0.  were 
added  to  each  sample  and  the  solutions  were  allowed  to  cool 
for  1/2  hour.   The  solutions  were  diluted  to  200  ml  with 
distilled  water.   A  titration  with  0.500N  FeSO.  was  done 
using  O-phenanthroline  indicator.   The  endpoint  is  reached 
when  the  dark  green  color  changed  to  maroon.   If  75%  of 
the  K2Cr07  was   reduced  by  the  soil  then  the  determination 
is  repeated  using  less  sample.   The  percent  organic  carbon 
is  found  by: 

M   K ?Cr07-M   FeSO,  x  0.3 

Organic  Carbon  =   eq  L '     -ec| *  x  1  33 

grams  of  dry  soil 


Problems  in  Analysis 

Many  problems  were  encountered  with  the  Hg  analysis. 
The  concentrations  in  the  soils  are  very  low  and  many  were 
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below  the  detection  limit.   The  standards  used  in  the 
analysis  were  not  stable  for  the  lower  concentrations  and 
had  to  be  prepared  every  2  hours  because  of  Hg  loss.   The 
solutions  were  somewhat  stabilized  with  HNO,. 

The  analysis  of  plant  material  was  difficult  with 
perchloric  digestions.   Hg  could  not  be  determined  in  sage 
samples  because  the  reaction  of  perchloric  acid  with  the 
finely  ground  sage  was  too  rapid  and  violent  resulting  in 
loss  of  Hg.   In  the  sage  determinations  it  is  suggested  that 
the  nitric  acid  be  allowed  to  oxidize  the  sample  overnight 
and  then  add  dilute  perchloric  acid  until  the  ratio  is  5:1. 
This  may  yield  better  results. 

Cadmium  was  also  determined  in  all  of  the  soil  samples 
but  the  results  seemed  too  high.   Some  random  analysis  for 
calcium  were  done  and  the  results  compared  with  Cd.   There 
was  a  high  correlation  and  the  Cd  results  are  probably  being 
influenced  by  the  known  Ca  interference.   Since  the  proper 
equipment  was  not  available  for  background  correction  the  Cd 
results  were  discarded.   Li  was  reliable  in  soil  analyses 
but  absorbance  could  not  be  distinguished  from  noise  in  all 
the  plant  samples. 

Boron  and  Molybdenum  Analyses 

B  and  Mo  analyses  were  done  at  Project  Central  Analytical 
Facility  (Meglen,  1977).   Soil  samples  were  analyzed  for  B 
by  the  following  method. 
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A  1  gram  sample  was  fused  with  2  grams  Na~CO,  in  a 
platinum  crucible  at  1000  C  for  1/2  hour.   This  was  cooled 
and  dissolved  in  HC1  and  diluted  to  50  mis.   All  labware 
used  was  polyethylene  to  avoid  contamination  from  B  in 
glassware.   A  1  ml  aliquot  is  buffered  with  ammonium  acetate 
EDTA  solution  to  a  pH  of  4.5.   The  analysis  is  done  colori- 
metrically  with  Azomethine-H  at  429  nm  on  a  Perkin  Elmer 
spectrophotometer.   The  analytical  error  is  5-101. 

Mo  in  soils  was  determined  by  placing  a  1  gram  sample 
in  a  100  ml  teflon  beaker  and  adding  concentrated  HNO,, 
H~S0.,  and  HF.   The  solution  is  evaporated  to  dryness  at 
240  C.   This  process  is  repeated  until  the  residue  is  white. 
One  ml  of  HC1  is  added  with  4  mis  of  deionized  H-O.   The 
solution  is  filtered  and  thiocyanate  solution  is  added.   The 
analysis  is  done  by  atomic  absorption  using  a  nitrous  oxide- 
acetylene  flame. 

The  B  and  Mo  in  plants  was  determined  by  ashing  3.00 
grams  at  450  C  for  24  hours,  and  carrying  out  analyses 
similar  to  the  soils,  except  that  Mo  was  determined  colori- 
metrically  in  plant  samples. 
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APPENDIX  II 

Geometric  Means  and  Deviations 

Geometric  means  are  calculated  using  standard  techniques 
on  log  (base  10)  transformations  of  data  except  for  pH  values 
which  are  already  log  values.   The  geometric  deviations  are 
analagous  to  standard  deviations  of  log  transformed  data. 

In  the  case  of  Hg,  part  of  the  data  are  below  detection 

limit  so  the  data  set  is  described  as  censored  (Miesch, 

1976a) .   The  detection  ratio  is  the  number  of  determinations 

above  the  detection  limit  (20  ppb  for  Hg)  divided  by  the 

total  number  of  samples.   For  Hg  the  detection  ratio  is 

222/253  =  0.88.   The  adjusted  mean  is  found  by: 

X  =  X1  -  X(X-X  ) 
v    oJ 

The  adjusted  deviation  is  found  by: 

s2  =  (s1)2  +  A(X'-Xo)2 

where  X  and  s  are  the  mean  and  deviation  of  the  uncensored 
data  and  X  is  the  detection  limit. 

Lambda  is  determined  graphically  (Cohen,  1959).   This 
method  reduces  the  uncertainty  in  the  calculation  of  mean  and 
deviation.   Dropping  the  less  than  samples  yields,  X  =  46.6 
ppb  and  s  =  1.56.   The  censored  calculation  values  are  X  = 
43.3  ppb  and  s  =  1.4. 

Further  adjustments  can  be  made  to  geometric  or  standard 
deviations  to  subtract  out  the  deviation  caused  by  analytical 
variability. 
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The  adjusted  geometric  deviation  (GDn)  is  calculated 


by: 


GD 


n  =  ^J(GD)2  -  (GE)2 


where  GE  is  the  geometric  error.   The  geometric  error  is 
found  by  calculating  the  geometric  deviations  of  all 
analytical  replications  (see  computer  program  Var  (Table  16 
in  Appendix  III).   The  geometric  error  is  the  average  of  all 
the  geometric  deviations  for  replications. 

The  use  of  geometric  deviations  is  most  important  in 
determining  95  percent  confidence  intervals  for  each  element, 

All  the  ranges  in  Table  2  have  been  calculated  using  GD 

°   n 

rather  than  GD.   This  narrows  the  range  and  gives  a  more 

realistic  estimate  of  element  concentrations. 

The  95%  expected  range  is  found  by: 

[ GM. j  tQ  GM  tQ  [GM  x  (    )1.96] 

(GDN)1'96  N 

where  1.96  is  commonly  rounded  to  2.0.   This  states  that  951 

of  the  samples  determined  by  this  study  or  any  other  in  this 

area,  should  lie  in  that  range  (assuming  analytical  methods 

used  are  similar  in  detection  limit  and  sensitivity). 


Analysis  of  Variance  Calculations 

The  computer  program  Nest  (Table  15  ,  Appendix  III)  was 
used  to  determine  the  estimated  components  of  variance  at 
each  level  of  the  nested  or  hierarchial  model.   The  problem 
in  the  interpretation  of  these  results  are  in  expressing 
the  validity  or  confidence  in  the  estimated  variance  compon- 
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ent .   Since  the  regional  component  is  the  most  important 
tests  of  the  significance  of  the  estimates  are  needed.   The 
standard  test  is  the  F-statistic  value.   As  can  be  seen  by 
inspection  of  results  very  few  of  the  highest  level  vari- 
ance components  are  significant  at  the  951  level.   This 
means  statistically  that  it  cannot  be  proven  that  the  com- 
ponent is  different  from  zero.   This  problem  arises  because 
of  the  distribution  of  the  degrees  of  freedom  in  the  nested 
model.   None  of  the  first  levels  (>1.6  km)  have  more  than 
three  degrees  of  freedom.   This  means  the  F-value  has  to  be 
very  large  if  the  variance  is  to  be  significant.   Since  the 
only  way  to  increase  the  degrees  of  freedom  at  the  highest 
level  is  to  collect  many  more  samples  and  hence  make  the 
cost  of  the  study  prohibitive.   The  values  of  variance  are 
accepted  as  the  best  estimates  available. 

Techniques  are  available  (Miesch,  1976b)  to  assess  the 
usefulness  of  variance  components.   The  variance  ratio 

XT  2 

N        s 

v  -  v  =      a 

D~    Z™  2X  Z~  2 
V   s,  +s  +s,  +  s 
b   c  d  e 

is  used  to  determine  the  efficiency  of  a  sample  design.   N 
is  the  variance  between  units  (sections)  and  D   is  the 
variance  within  units.   With  the  variance  ratio,  the  effec- 
tive number  of  samples  collected  at  random  (N  )  can  be 
determined  graphically  (Miesh,  1976b).   With  N  ,  the  maximum 
permissible  error  variance  of  the  means  can  be  calculated: 


T  2009  59 


e2x  2^  2A  2 

S,  +5  +  S,  +  S 

r,   _   b   c   d   e 

E   -  n 

r        Nr 

The  E   refers  to  a  balanced  design,  the  maximum  error  vari- 
ance for  a  nested  design  is  given  by: 

2     2         2  2 

S,         S  St  S 

E   =  J.   +  „  c„   +  „ d +      e 


s    Nb    Nb-Nc    Nb.Nc.Nd   Nb.VNd.Ne 

where  N,  is  the  number  of  randomly  sampled  areas  within  each 

section  and  so  forth  for  N  ,  N-,,  N  . 

c'   d'   e 

Another  useful  quantity  is  the  variance  mean  ratio  (V  ) , 

V   =  N  /E 
m    v'  s 

for  a  non-hierarchial  design  E   is  used  in  place  of  E  . 

b        r  Y  s 

At  this  point  an  example  of  these  calculations  is  neces- 
sary.  For  molybdenum  in  soils  the  computations  are: 

s2  =  0.01485  s2  =  0.0159 

a  d 

si    =  0.0  s2  =  0.00273 

b  e 

s2  =  0.04625 
c 

The  variance  ratio  (V)  is, 

v    0  +  0.04685  +  0.0159  +  0.00273    U'^ 

from  graphs  N   =4 

-     0  +  0.04625  +  0.0159  +  0.00273 

E   =  -. 

r  4 

E   is  found  by: 


Es 4 ('91)  +  4T8  +  4-8-16  (-98)  +  4.8-16.32 


0.00273 
4- 8-16-32- 6 
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E   =  0.0035 
s 


where  0.91  and  0.98  are  correction  factors  explained  in 
Miesh  (1976a).   The  variance  mean  ratio  (V  )  is: 


V. 


_.  0.01485 
m  "'  0.0035 


V   =  4.28 
m 


From  experience  and  computer  simulation  studies  it  has 
been  shown  (Miesh,  1976b)  that  if  E   is  greater  than  E   then 
the  model  produces  at  80%  confidence  the  variance  between 
the  studied  units  (sections) .   The  801  is  derived  from  the 
graphs  of  N  .   In  other  words  if  the  maximum  error  of  the 
nested  design  (E  )  is  less  than  the  maximum  error  variance 
of  a  balanced  design  (E  )  then  the  model  chosen  signifi- 
cantly displays  the  differences  in  the  study  area.   These 
same  calculations  can  also  be  used  to  determine  the  most 
efficient  sampling  design  to  describe  variance  not  described 
by  a  pilot  or  initial  sampling  study.   If  it  was  found  that 
E  >E   for  many  components,  then  the  E   expression  could  be 
altered  until  the  observed  error  variance  was  acceptable. 

The  variance  mean  ratio  (V  )  serves  as  a  measure  of  the 

m 

stability  or  reproducibility  of  a  geochemical  map.   If  V 

m 

=  1  basic  differences  can  be  shown,  if  V  >3  then  the  result- 

'     m 

ing  map  should  be  quite  stable  (Miesh,  1976b).   As  can  be 

seen  in  Table  2,  all  the  V  values  are  greater  than  or  equal 

to  2.9  and  almost  all  the  E  values  are  less  than  E  .   This 

s  r 
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concludes  that  the  sample  design  was  sufficient  to  describe 
the  variance  between  units  and  that  the  maps  should  be  very 
stable. 


Student' s  t  Test 

To  determine  the  significance  of  differences  between 
two  means  the  Student's  t  test  is  used: 


t  = 


Xl  "  X2  „ 

2 2"' 

(Sj)     (s2) 


N. 


N. 


C5lX  +  s22/N2)2 
(s1)2/N1    (s2)Z/N2 
M,-l —  +   N„-l_ 


'1      "2  -1  2 

This  does  not  account  for  differences  in  variances  for 
each  population  tested,  which  is  not  of  great  concern  with 
the  data  in  this  study.   Student's  t  test  data  show  high 
significance  in  many  of  the  populations  studied. 

Other  Data  Reduction  Techniques 

Linear  correlation  coefficients  are  calculated  for  the 
log  (10)  transformed  data.   Log  transformation  is  necessary 
because  the  parts  per  million  data  would  give  spurious 
correlations  due  to  excessively  high  and  low  values  in  the 
population. 

Multiple  regression  analyses  were  also  run  on  the  data 
to  determine  how  each  element  varies  with  the  other  para- 
meters.  Trend  surface  analyses  were  also  done  to  display 
regional  and  lithographic  influences  over  the  grid  area. 
The  correlation  coefficients,  regression  analyses,  and  trend 
surface  analyses  results  are  supportive  of  each  other  in  the 
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interpretation  of  the  data.  All  three  of  the  above  methods 
were  performed  on  library  computer  programs  at  the  Colorado 
School  of  Mines . 
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Table  11.   Analysis  of  Variance  Results  of  Soil  Samples 


Source  of 
Variation 


Sum  of 
Squares 


Hg 

1.6  km 
0.4-1.6 
0.1-0.4 
0-50  m 


km 
km 


0865 
2127 
4949 
152 


Degrees  of 
Freedom 


3 

4 

8 

16 


Mean 
Square 


F 

Value 


Estimated  Percent 
Variance   of  Total 
Component  Variance 


Replications  0.2121 

Zn 

1.6  km      0.0504 


0 
0 
0 
0. 


0288 
0532 
0618 
0719 


0 
0 
0 
3 


54 
86 
86 
39* 


10 


0.0212 


0.0168   1.82 


0 
0 
0 
0 


042 


0.021 


0.0007 


0 

0 

0 
66.5 
33.5 


7.6 


0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

Li 

1.  6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

B 

1.6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

0.0369 
0.1263 
0.0719 
0.0149 

4 

8 

16 

8 

0.0092 
0.0158 
0.0047 
0.0019 

0.584 
3.36* 
2.53 

0 

0.0045 
0.0024 
0.0019 

0 

47.4 
25.3 
19.6 

0.6232 
0.2505 
0.6077 
0.0721 
0.0198 

3 

4 

8 

16 

12 

0.2077 
0.0626 
0.0759 
0.0045 
0.0016 

3.31 

0.82 

16.85* 

2.73* 

0.0130 

0 

0.0278 

0.0023 

0.0017 

29.0 

0 
62.2 

5.1 

3.7 

0.0258 
0.0965 
0.5425 
0.1379 
0.0191 

3 

4 

8 

16 

11 

0.0086 
0.0241 
0.0678 
0.0086 
0.0017 

0.36 
0.36 
7.86* 
4.96* 

0 

0 

0.0221 

0.0056 

0.0017 

0 

0 
75.2 
18.9 

5.9 

Mo 

1.  6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

As 

1.  6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

0.4606 

0.0483 

1.016 

0.3351 

0.0164 

3 
4 
8 
16 
6 

0.1535 
0.0121 
0.1270 
0.0209 
0.0027 

12.7* 
0.09 
6.06* 
7.66* 

0.0148 

0 

0.0463 

0.0159 

0.0027 

18.6 

0 
58.0 
19.9 

3.4 

0.7521 
0.0307 
0.2534 
0.1035 
0.0848 

3 

4 

8 

16 

19 

0.2507 
0.0077 
0.0317 
0.0065 
0.0045 

32.64* 
0.24 
4.89* 
1.45 

0.0193 

0 

0.0083 

0.0014 

0.0045 

57.7 

0 
24.9 

4.0 
13.4 

Organic  Carbon 

1.6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

1.469 

0.2886 
0.4549 
0.5454 
0.0275 

3 

4 

8 

16 

10 

0.4895 
0.0722 
0.0569 
0.0340 
0.0028 

6.78* 
1.27 
1.67 
12.39* 

0.0398 
0.0024 
0.0077 
0.0255 
0.0028 

50.9 
3.1 
9.9 

32.6 
3.5 
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Table  11.  (Continued) 

Estimated   Percent 
Source  of   Sum  of   Degrees  of   Mean     F    Variance   of  Total 
Variation   Squares    Freedom    Square  Value   Component   Variance 

EH 

1.6  km 
0.4-1.6 
0.1-0.4 
0-50  m 


0 
1 
2 
1 
Replications  0 


km 
km 


0974 

048 

448 

644 

2267 


3 

4 

8 

16 

20 


0329 
2619 
3059 
102  8 
0113 


0 
0 
2 
9 


13 
86 
98* 
08* 


0 
0 
0 

0 

0 


0621 
0597 
0113 


0 

0 

46 

44 


Analytical  Variability  Determined  from  All 
Samples  Analyzed. 


Element 

Hg 

Zn 

Li 
B 

Mo 
Org  C 

pH 


Percent  Deviation 
17.7 
5.9 
3.1 
5.1 
8.6 
5.8 
0.8 


D 
0 
D 
G 
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Table  12.   Analysis  of  Variance  Results  from 
Sagebrush  Samp 1 e s . 


Estimated 

Percent 

Source  of  Sum  of 

Degrees  of 

Mean 

F 

Variance 

of  Total 

Variation  Squares 

Freedom 

Square 

Value 

Component 

Variance 

Zn 

1.  6  km 

0.4306 

3 

0.1514 

1.31 

0.0037 

2.1 

0.4-1.6  km 

0.4624 

4 

0.1156 

1.67 

0.0103 

6.0 

0.1-0.4  km 

0.5527 

8 

0.0690 

0.42 

0 

0 

0-50  m 

2.586 

16 

0.1616 

1.23 

0.0265 

15.  5 

Replications 

B 

1.  6  km 

0.9182 

7 

0.1312 

0.1312 

76.4 

0.0861 

3 

0.0287 

2.87 

0.0022 

31.6 

0.4-1.6  km 

0.0399 

4 

0.0099 

2.06 

0.0011 

16.5 

0.1-0.4  km 

0.0387 

8 

0.0048 

1.80 

0.0010 

14.9 

0-50  m 

0.0403 

15 

0.0097 

11.07* 

0.0023 

33.5 

Replications 

0.0007 

3 

0.0002 

0.0002 

3.5 

Mo 

1.6  km 

0.4306 

3 

0.1435 

3.82 

0.0100 

18.3 

0.4-1.6  km 

0.1504 

4 

0.0376 

0.39 

0 

0 

0.1-0.4  km 

0.7698 

8 

0.0962 

4.39* 

0.0276 

50.2 

0-50  m  ' 

0.3283 

15 

0.0219 

8.04* 

0.0146 

26.6 

Replications 

0.0354 

13 

0.0027 

0.0027 

4.9 

'Significantly  different  from  zero  at 


a  ■ 


05 


Analytical  Variability  Determined  from  All 
Sage  Samples  Analyzed. 


Element 
Mo  (sage) 


Percent  Deviation 
4.9 


B  (sage) 


2.0 
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Table  13.   Analysis  of  Variance  Results  for 
Ricegrass  Samples. 


Source  of 
Variation 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 
Square 

F 

Value 

Estimated 

Variance 

Component 

Percent 
of  Total 
Variance 

Hi 

1.6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

Zn 

1.6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

B 

1.6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

0.1618 
0.1502 
0.1395 
0.4051 
0.0268 

3 

4 

8 

13 

4 

0.0539 
0.0376 
0.0174 
0.0312 
0.0067 

1.44 
2.15 
0.56 
4.65* 

0.0018 

0.0049 

0 

0.0222 

0.0067 

5.1 
13.8 

0 
62.3 
18.8 

0.1454 

0.3047 

0.6943 

3.386 

0.0906 

3 
4 
8 
13 
3 

0.0485 
0.0762 
0.0868 
0.2605 
0.0302 

0.64 
0.88 
0.33 
8.62* 

0 

0 

0 

0.2138 

0.0302 

0 

0 

0 
87.6 
12.4 

0.0829 
0.0918 
0.2206 
0.1735 
0.0056 

3 
4 
8 
16 
7 

0.0277 
0.0229 
0.0276 
0.0108 
0.0008 

1.21 

0.83 

2.54 

13.63* 

0.0004 

0 

0.0068 

0.0089 

0.0008 

2.1 

0 
40.4 
52.8 

4.7 

Mo 

1.6  km 
0.4-1.6  km 
0.1-0.4  km 
0-50  m 
Replications 

0.2832 
0.0915 
0.3785 
0.2972 
0.0587 

3 

4 

8 

15 

10 

0.0944 
0.0229 
0.0473 
0.0198 
0.0059 

4.13 
0.48 
2.39 
3.38* 

0.0070 

0 

0.0099 

0.0117 

0.0059 

20.2 

0 
28.8 
33.9 
16.9 
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Table  14.   Analysis  of  Variance  Results  for 

Wheatgrass  Samples. 


Source  of  Sum  of   Degrees  of 
Variation   Squares    Freedom 

Mi 

1.6  km 
0.4-1.6 
0.1-0.4 
0-50  m 


Mean 
Square 


Estimated 
F    Variance 
Value   Component 


km 
km 


3 
4 
8 
13 
4 


0.1618 
0.1502 
0.1395 
0.4051 
Replications  0.0268 

Zn 

1.6  km       0.0206  2 

0.4-1.6  km   0.0343  3 

0.1-0.4  km   0.0981  6 

0-50  m       0.2945  11 

Replications  0.0859  5 

B 

1.6  km  0.0126  2 

0.4-1.6  km  0.0201  3 

0.1-0.4  km  0.1872  6 

0-50  m  0.2443  10 

Replications  0.0054  5 

Mo 

1.6  km  0.0484  2 

0.4-1.6  km  0.1416  3 

0.1-0.4  km  0.3529  6 

0-50  m  0.0946  10 

Replications  0.0483  5 


0 
0 
0 
0 
0 


0 
0 
0 
0 
0 


0 
0 
0, 

0. 
0, 


0 
0 
0 
0 
0 


0539 
0376 
0174 
0316 
0067 


0103 
0142 
0163 
0268 
0172 


0063 
0067 
0312 
0244 
0011 


0242 
0473 
0588 
0095 
0097 


1.44 
2.15 
0.55 
4.65* 


0.90 
0.69 
0.61 
1.56 


0.94 

0.21 

1.27 

22.42* 


0.51 
0.80 

6.22* 
0.98 


0 

0 
0 
0 
0 


0 
0 
0 
0 
0 


0018 
0049 

022  2 
0067 


0 

0 

0 

0.0088 

0.0172 


0023 
0203 
0011 


0 

0 

0.0230 

0 

0.0096 


Percent 
of  Total 
Variance 


5.1 
13.8 

0 
62.3 
18.8 


0 

0 

0 
33.8 
66.2 


0 

0 

9.8 

85.6 

4.6 


0 

0 
70.4 

0 
29.6 


H 
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APPENDIX  III 

Computer  Programs 

The  program  Nest  (Table  15)  was  used  to  examine  the 
analysis  of  variance  samples.   Figure  22  shows  the  hier- 
archial  model  with  three  input  parameters  (NL's)  that 
describe  the  nested  sampling  model  used  in  this  study. 

The  programs  used  to  find  geometric  means  and  devia- 
tions (Table  18)  and  analytical  variability  (Table  16)  were 
adapted  from  Klusman  (1976) . 

Data  Files 

Table  19  lists  the  analysis  of  variance  data  and  the 
plot  is  shown  in  Figure   6.   Table  21  lists  the  grid  data 
and  Figure  5   shows  the  plot  of  this  data.   The  flags  used 
in  both  data  sets  are  on  the  line  proceeding  the  data  for 
each  point.   A  zero  indicates  the  sample  was  collected  on 
the  Uinta  Formation.   A  90  shows  the  sample  was  taken  on  the 
Parachute  Creek  member.   A  08  shows  that  the  sage  sample 
was  subspecies  wyomingensis . 

Figures  18  to  21   show  the  data  used  in  hypothesis 
tests  (Groups  1-4). 


H3 

IV) 

o 

O 


1 


4     A 


A    4 


*   a 


Figure22.  Analysis  of  Variance 
Model  Showing  Values  of  NL  used 
in  the  Computer  Program  Nest. 
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1  woo'    "ssav^sls  &r- used " caloulate     70 

iC   ""oftSJft  Eft!""!."."  ",Si:  N0r,M.,.C!««9.t..I««» 

1   JlipNSlON  N<?<999).YV(999>,YS<999) 

0  PENSION  A^A(999).SIC(999).PCSICC999.) 
Ol^ENSIOV  WMSi(999),rvAL(999) 
COUBLE  PRECISION  A.3iC 

J      rlllWnV.'V^   !^uT  file  N»mE  for  nls 0»fl 

READ  (4.6)  A 
5         rORMAT(lA!0) 

,       ^."iSl-ENTER  INPUT  FILE  NA«  FOR  0*TA  VALUES OaT-) 

rEaQ(4,6>  B 

8       fOR*5t?15! 'ENTER  OUTPUT  FILE  NAME  I  DAT') 

REaO(4.6>  C 
§PEN(UNlT=3,riLE»C) 

OPEN<UN'lT  =  l-.f  ILE  =A> 

0PEN(UMT  =  l2triLE  =  3) 

IFLAG=Z 
NCT  =  1 

nYCTs0  „    __ 

IF    (IFLAG.E2.35G3    TO    52 
J.Z  IF(IFLAG,EQ.»)    GO    TO    4 

5        SSJmt      N    F    R      XAM^I,     '    E "a    Ful'oF    99    AT    T^END    OF    EACH 
C  UMr\TA TI3 N    F CK       JAM    u-  Tw  nL    SETf     JF    Ng 

J        -rJtNG  OF    T    r    n[   W-MTS    II   OESIRCO    A    FLAG    0F    93    I«    SUFFICIENT. 
\        T<?M   EHTMcSSJ   PC"    tINE    SEPARATED    BY    -    ,    -    BEGIN    A    Nd 

C        LINE    AFTER    a    93    CR    99    FLAG   HAS    BEEN   ENTERED. 
2-      -  REA:<1Z,32)(^-L(n'l=l'2?) 

3;  rORMAT(2ZI) 

qO    40    1*1.22 

IF    <NL(I) .EQ.90)    GO    TO    46? 

NS{NCT)«*L{I> 

NCTsNCT*l 

4:  CONTINUE 

C        TYPE^N^fvALUE    OF    NF  I    NUHBER    OF    FACTORS    IN    THE    MODEL 

'r  ^R^a^.'EUTER    THE    VALUE    FOP    Kf/I 

READ(4,7?)     NF 

7-      ro^MAT(U) 

WF1=\F*1 

Nr2  =  -.F»(\F*D/2 

c   ».ea-jNin"w  first  SET  Of  NL'S 
m-    pead(i7,9:)(hl(I>»  m.2c> 

9;       rORMAT{?OI) 

no  122  I«l. 20      T„ 
ir(fl'L{i).NE.98J  GO  TO  103 
iFLASsi 
GO  TO  12^ 
i:/      jr(Wu( l) .HE. 99)  GO  T0  113 
GO  TO  125 
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11J     nSCICTjs  NL<!) 

\'CTs'jCT*1 
12.;     CONTINUE 

GO  TO  82 
C    hCil  IN  THE  rATA(Y),  (OME  PIECE  PER  LINE) 

c    usr  4  flag  or  Q9  *t  r*z   E-o  or  the  set 

C    CO'VERSIVJ  Of  OAT;  TO  L3G(18)  VALUES  IF  OESIREQ 

125  -RITE(4,126)  ,.  .  _ 

126  rOR'UTdX,' ;r  LOG  CONVERSION  Or  OATA  IS  OESlREO  TYFE  2i  IF  NO' 
1   TYPE  iv) 

PEA^U,??)  ilcg 
13.*     .  Jal 
14:-      REAT<12»152)  YY(J^ 

lF(YV(j$ .E0.99.)  CO  TO  19? 

irc!L0G.E3.D  GO  TO  145 

Y<J)  =  AL0G1Z  (YY(J)) 

GO  T0  162 
145       V(JJ»(YY(J>) 

152    .ro3HAT(ri3.3) 

164      KYCTsNYCT*l 

YSCiYCT)!  Y(J) 

J«J*1 

GO  TO  147 
C    -RITE  OUT  C*TA  VALUES  TO  TTY  AND  FILE 
182      WRITE  (8.135) 
185      FORHAT(lX,//ii:X, 'DATA  AND  NL  VALUES'.//) 

WRlTE(5.<82)<Y(J)rJ=l.NYCT) 

WRIT£(8,233)<NS<  I>.I=1.NCT) 
WRITE(4,192)(NS(D.I31»NCT) 

193  FORv'ATf  /22I2/) 
222  F0RMAT(/*3!2/) 
C    CALL  ANTSTu  SUBROUTINE  TO  ANALYSE  NESTED  0*TA 

CALL  ANESTy(\r,NS.Y,GM.S,NOF,EMS. I  UK,  IER) 
C    CALCULATE  k£An"S£uaRE  VALUES 
•  00212  1=1. Nf 

WHSO(I)=  SU)  /  FLOAT(NDFU)) 

212      CONTINUE 

C     CALCULATE  r  STAT,  VALUES 

DO  220  I=1.NF-1 

FVALCI)  •  wmsq(I)  /  WMSQ(I*l) 
222     CONTINUE 

wRlTru.233) 

WRITE(3,23?) 
23,'  F0RVAT<1*, //'EJECTED    MEAN    SOL»RE    COEFICJ.EMTS  •// > 

24,;  F0R,,aT(F1.t.5.  115) 

IS  =  1 

irs-;F 

C    *"?ITE  0;iT  EXPECTED  >"EAN  SQUARE  COEFFICIENTS 
HO  25?  N»l#vF 

w.RlTE(4,263)  <EUS<!  )  .  I  =  IS,  IF) 
v,RITE(s,26Z)  (E-'St  I  )  .  Is  IS,  IF) 

:s=if»i 

!FsIS*(\F-M)-l 

252     CONTINUE 

26J      F0R"ATfi2ri2,4) 

*.RITE(4,27?) 
C    CONFIRM  VALIDITY  QF  JATA,  IF  BAD  PRGm,  STARTS  OVER. 
27"?      FOfl,,AT<iX,»!F  -ATA  IS  GOOD.ENTERl,  IF  BAD.EnTER  3'  >_ 

hEa'»  (4,292)  JPRINT 


T2009 


72 


293 

m 

y 

2' ; 

31.' 

32  J 
33' 
C 
343 


353 

C 

363 

373 


e 
c 


38: 
C 


39J 

C 


«1<? 

*23 


i 
r, 

type: 

r 

R 
F 

w 
F 
F 

CALCU 

s 
s 
1 
s 

1  ) 
I 
s 

1    - 

! 

s 
1 
I 
I 
s 
1  - 
1    ) 

F 
WRITE 


1 

1 


o  n  i 

RjTEU 

IN  A  NO 
OR'UT( 

EADU, 

OR'UT( 
RITE  ( 

0fl»«4Tc 

ORH*T( 

LATE  c 
IG£MF) 
IG(MF- 
F  C«NF 
IGJNF- 
>/F.MS< 
FC(NF- 
IG(\r- 
D-EHS 

f<cnf- 

IG(\F- 
)-CMS( 
0) 

F((NF- 

IGfNF- 

D-EMS 

-EMS(N 

0R>'4T( 

OUT  C 


NT.EQ.D  GO  TO  29? 

Z 

,3J^) 


wR 
1*. 

31C 

4A5 

a. 3 

1*2 

1«, 
STI 
=  S< 

1)  = 
-2) 


!TE 

I  EN 
)  (A 
) 

23) 

//, 

MAT 

\r) 

<S< 

.EC 


2>={S( 
NF2-3) 
3J.EQ. 
3>=C3< 

(NF2-7 
4).EQ. 

4)=<S< 
NT2-12 

5J.ET. 

5)=CS< 

CNF2-1 
F2-163 

1X.F15 
OLUMN 


OUT  AREA  TITLE 
TER  AREA  TITLE, UP  TO  24  SPACES'/) 
REA(L).L=1,4> 

(aREA(L)  ,L=1,4) 
A5>/> 

3X,,T0UL  SUm  OF  SQuARES'.SX. 'DEC.  OF  FREEDOM/) 

E?  COMPONENTS  OF  VARIANCE 

/"OF(NF) 

\F-1)/Nl)F(NF-1)  -  SIG(NF))/EMS(NF2-1) 

,2)    GO  TO  363 

NF-2)/N0F<NF-2)-SlG(NF)*EMS(NF2-4)«SlGCNF*l 

3)G0  TO  362 

NF-J)/  ^0F(»jF-3)-Slc(NF)-EMS(NF2-8)«  5IG(NF 

).SIG(NF-2))/  EMS(NF2-6) 

C)  GO  TO  360 

Nr-4)/N0F(NF-4)-SlG(NF)-EMS(NF2-l3)*SlG(KF* 
)#SIG(NF-2)-EM3(NF2-ll)*SIG(NF-3))/EKS(NF2- 

2)    GO  TQ  363 

k.F-5)/NDF(NF-5)-SlG(NF)aEMS(NF2-19)«SIG(NF- 

S)»3!GCUF-2)-EMS{!NF2'l7)*SlGl.NF-3 

•  SIGCNF-.4J  )/  EMS(NF2-15) 

,6,2Xi  U8«FiC,6  #5XtFlPl,6   //) 

HEADINGS  FOR  ANOVa  TABLE 


WRlT£(a,373) 

WRITE(4037?) 

FORKAT(iX, 'LEVEL' >10X, 'SUM  OF  SQUARES  ', leX ,• OEG .  CF  FR,», 
17X,'MEAN  SCUARE'.lCX.'F  VALUE', ax, 'EST,  COmP,  OF  VAR.'„ 
1   6X,'PCT.Cf  TCT.VAR.'///) 
CALCULATE  TOTAL  ESTIMATED  COMP,  OF  VaR  I  ANCE  (  tOts  I  G  )  OMMTING 
VALUES  LESS  TM4N  qr  EQUAL  TO  ZERO 

D0333  1=1, NF 

IF((SIGCD)  .LEV  3.)  GO  TO  383 

TOTSIG=  TQTSIG*SIG<!> 

CONTINUE 
CAlCULAtE  PERCENT  cF  total  VARIANCE  (PCSIG) 

00393  1=1, ^F 

PCSIG(T)»  (StG(l)  /  TOTSIG  )  •<1?3) 

!F(fSiG(  I))  iLE.  C)  PCSIG(  I)=3, 

CONTINUE 

F0R"UT(F1?,5) 

--Rite  cut  ancva  ta^le 

D042J  1=1, nf 

wRPE(8t*l?)  I,S(D,NDF(I),WMSQ(I),FVAL(I),SIG(I).PCSIG 
KI) 

.OlTEU,413)IfS(I).N0r{I),WMS0(I),FVAL(I),SlG(I),PCSIC 
Id) 

FOR"AT(iX,3XiI3illX.Fie.6#l3X,I3.13X,Fi0,6,8X,FlZ.4.12X 
1   ,Fl.T.5,lCX,Fia.4i///) 

CONTINUE 

WRIT£C4,332> 

WRITE(8,332) 

-RlTE(4,A4T)S(fri) # NOF(MFl) 

#.RlTf    TOTAL    SUM    Or    SQUARES    A  JO    DEGREES    OF    FHEEQQM 
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c   -rite  grand  mea*  a*-o  error  parameter 

wRITE(4,43Z) 
wRlTF<8,432> 

43/      rORMAT(ix,//,2^.  '  CRA  jD  MEA'J'tlBX, 'ERROR  PARAMETER'/) 

WRITE(8,«4?)  G^.IER 

*RlTEC4,44?)3y,  US 
44.'      F0R"4AT(riZ,3,l-5X,!l":) 
c   E'^Ea  e  I"  MCRr  :ata.  1  TO  Mr.RGE  OaTa  SETS,  2  TO  STCP 

URITE(4,452) 
45,}     rORMATciX,  'If  "ORE  DATA, ENTER  2,lF  MERCE  ENTER  It  IF  STC 

IP  ,  ENTER  2') 

REA?»<4,2P0)  iflag 

IF(IFLAG.EQ.2)  CALL  EXIT 

GO  Tfl  i? 
C    ANALYSE  MERGED  DATA  SETS 

C    mDTE  A  "EW  SE?  OF  \L'S  MUST  3E  ENTERED  CONSTRUCTED  FRCf 
C    THE  NEWL*  COMBINED  OR  MERGED  VALUES 
463      WRITE(4.63) 

READ«4.7?)  NF 
470      FORHATc (131)) 

wRITE<8.472>  (VS(PiI«t*NCT) 

WRITE(8.<8Z)  (VS(t). J*1»NYCT) 
483      F0RMAT(8F12.4) 

NFl«NF»l 

NF2«NF«{NF*i)/2 

CALL  ANESTU(^F,\S.YS|GM,S,NDF,EMSt IWK,IER) 
C   ENTER  ELEMENT  NA"E 

WRITE<4,492) 
490     F0R'-UT<iX,  'ENTER  ELEMENT'/) 

READM,5?2)  ELEM 
533      F0R"UT(A6) 

URlTE(e,5lZ)ELEM 
510      F0Rl'AT(iX,A6) 

WRITEC4,24Z)(S(IJ»N0F(I),!=1,NF1) 

wRlTEU.262) (EHSt I) . 1=1. NF2) 

*RITE(8,24?)  (S(H»fOF(I),I=l,NFl) 

wRITE(B,26^) IE-S(I) . t=l.NF2) 

WRITEC8.443)  G^.IES 

IFLAG=2 
GO  TC  347 

STCP 

END 


[ 


fA. 

Table  16.  listing  of  Program  fgr  Analytical 

Variability.      .        ._  _    .  .  .- 
C  P»5G«4M  TO  DETERMINE  A'lALVTfCAL  VARIABILITY  OF  REPLICANTS 
C  ^ATA  IS  RE*0  rSQM  Fqrjs.OaT  ,  FORmaT  MyST  9E  Entered  fcr  each  element 
C  IH  STATEHT^T  *  173.  THIS  VERSION  WILL  ACCOMODATE  DATA  WHERE  SCME  CR 
C  ALL  SAMPLES  ARE  REANALYZES!  ALSO  IF  SOME  ELEmEnTs  ARE  NCT  QETER^NEC 
C  F3-'  ALL  SAMPLES. 

DIMENSION  OATA(lT) 
b    aEAO(19,ie?)  LABEL 
l?j  rnf,MAr(A<) 

IF(LA3EL.E2.4W9999)  GO  TO  1 

SSUM1»J, 

GSuM2  =  "<. 

GSUM3=2. 

<»£ 

^ITE(6.131J  LABEL 

101  F0RMAT(lHli32X,A4) 
wRlTE<6#l?2> 

102  F0RMAT(iH2,4HQAT4/56H  ARITHMETIC     STANDARD    PERCENT   GECMETRIC 
1   GEO^ETRIC/54*       MEAN     DEVIATION   DEVIATION    KEAN     DEvUT 

2I3N) 

7  CONTINUE 

DtJ  11  1=1  #3 

READ(8.1?3)CATA( I) 
123  FORMAT(54X,f5.1) 
11       CONTINUE 
N»l 

SUm«0. 

s'JM2=e. 

SUHL«?. 

C  IF  ThE  FIRST  VALUE  OF  DATA(I)  =  0.3  REaD  THE  NEXT  VALLE  SINCE  THIS 
C  ELEMENT  WAS  NOT  DETERMINED  FOR  THIS  SAMPLE 

IF(DATACI)  .*En.?.>GO  TO  7 
4  IF<DATA(N) ,EQ.x\ )  GO  TO  2 

IF(DATA{\)'.EC  9.9)  GO  TO  3 

lF(DATA(N).LT,3,r?3)  DATA < N )= ASS < DATA < N ) > 

SUM  =  SUM*daTa('J) 

SUML  =  S'JML*ALOGlZ(DATA(N)  ) 

SIJm2  =  SJM2*(DATA('i)  )*»2 

SUML2  =  3UML2*(AL0G1-'<0ATA(N)  )»»2) 

VaN  +  1 

50  TO  4 

2  M*N-1 

-:»lTE(7.ie4)(DATA(I).I=l»N) 

134  F0RMATC1H  il?F6.?J 
X'lsFLCA-  N) 
XMpAN=aUM/XN 
GMFANSFXP{2.3J3»(SUML/XN>> 

iFCi.EC.DGO    Tn    99 
XDEV  =  S:RT((1./((xN»*2)*'XN))»((  <X\»SUm2)-«SuM»»2)  ) ) ) 
j-5EV  =  EXP(2.333»(SQRT(  < 1 ./(( XN»»2 ) -XN >)•<{< XN»SUML2 >- C SUfL«»2 >))>> ) 
GO    n    fl5 
99  X0EV»2 

C!)EV«3 
85  SPE=D=0, 

0*>    5    1  =  1.  K 
5    5PER0  =  SPERD*ABS( ( (DATA ( I J-XMEAN) /XmEAN)  »10e  .  ) 

-EnO=S?E«D/XN 

^RlTE(6.i:5)  X-EaN,XOEV,PERD,GMEaN,CDEV 

135  f3RMAT(lH    »5FU.3> 

C  IF  DaTa<2>*0.3  IT  IS  ASSUMED  ThaT  THIS  SAMPLE  WAS  NQT-RE- — 


I 
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C   ANALvIfO    a\!0    TU£    sirxf    SAMPLE    SHOULD    BE    READ    SO    Twe    COLDER    V.  ILL    N^T 

C  =5E  altereh  ano  the:  correct  avg.  deviations  can  be  computed. 

IF(6aTa(2).EC.2. )CO  TO  7 

GSuMl«GSUfl*XDEV 

CSUH2  =  0SUW2  +  PE^'0 

GSuM3=GSUM3*GDEv 

K«K*1 

GO  JO  7 
3  XKsfLOAT(K) 

AVOEVsCSU^l/XK 

-RITE(6,1C6)  LABEL. AVOEV 
126  rORMAT(iHP,3l^AVERAGE  STANDARD  DEVIATION  FqR  .A4.3M  =  ,Fte,3) 

*VPER  =  GSUH2/X'< 

wPITE(6i127)  LABELiAVPER 
107  FORHATdw  .33HAVERAGE  PERCENT  DEVIATION  FOR  ,A4,3>"  *  ,Fl3,3) 

AGDEV«GSUM3/Xi< 

wRlTE(6,ie5)  LABEL. AGDEV 
138  FORMATCiH  ,32^AVERAGE  GEOMETRIC  DEVIATION  FOR  ,A4#3H  s  ,FlC,3> 

GO  TO  6 
1  END 


T2009       Table  17  •  Listing  of  Program  used. to  76 

Perform  %  Teats. 

'TENSION  LAPEL (12)  ,  irMT(i6)l!HEA0(t6)fJHEA0(l6),CAT*l(999ll2)f 
1?ATA2C999»1?).  1QPT(12)iX2AT1(999)  ,xO*T2(999)  ,ICRCE:K(323)  iISa^PK 

1  3c3>»ISamP2(3;:d 

C  READ  IN  NU*gFR  or  DaTA  BATCHES,  NUMBER  Or  ELEMENTS 

C  aBaTo  is  DEEPED  as  2  SETS  Or  DaTa  FOR  WHICH  THE  T-TEST  FOP 

C  2lFFEOE»iCj;  IN  H£aNS  IS  BEING  RUN 

PEaO(13.1?7)  uy,uMM 

IS?  rCRMAT{2I2) 
C  *£AD  IN  ELEMENTS  DETERMINED 

=>EaD(13,101) ( LABEL ( I) . 1=1 »12 ) 

101  FTRMAT(12A4> 
C  READ  In  LD3  OR  NORMAL  OPTION  FQR  EACH  ELEMENT 

REA0(13.lZ8)(I0PT(I)#I=lil2) 
108  F0RMAT(12A3) 
C  3£aD  In  PLaG  FOR  END  OP  D*TA  SET 

PEaDC13»1Z21)  ^LAG 
1721     F0RMAT(1?F3.?) 
C  3EJ5IN  LOOP  THRC'jCH  EACH  BATCH  OP  DATA 

DO  1  Mil.M- 
C  *SaO  In  HEADING  FOP  FACH  DATA  SET 

BEAD C 13.133) (IhEaD(I). 1=1.16) 

REaOC 13,123) (J^eaOCD » 1=1.16) 

31  FORHAT(l3,l7XfF4,3,lX,P3.e,lX.F3.0ilX.F4r2,iX,F4,l,lX,F4,2,lXi 

1   F5.:.1X.P5,1.2X,F4,2/) 

9999  FORMATdX,  I  3  .  9  (  3X  ,  F6  ,  2  )  ) 

36  FORUT(n.l7X,F4,0.lX,F3.0,lX.F3.0.lX.P4,2,iX,F4,l.lX,P4,2, 

1    1X.P5.3.1X,F5,1.2X.P4,2/) 

Nl«l 

N2  =  l 

32  F0RV'ATM2) 

3  REaDO.32)  ICeEEi<tI) 
IP(ICREEK( I) .E2.93)    GO    Tq    3 
IP(ICREEK(I),E0.9)    GO    TO    993 
ir(ICRFEX(I).E3.3)    GO    To    3 

READC8.31)     ISa'-'P1CN1)»(0aTa1(N1,K),K  =  1,MmM) 

WRITEU.9999)  ISAMP1<.'«1>  ,  (OATAHNl.K)  ,K*1,HMM) 

ir(0ATAl(Nl.l).E3.PLAG)  GO  TO  2 

M  =  Nl*i 

GO  TO  3 
2       M^i-1 

GO  TD  199 
993      REA'(S.36)ISAvP2(N2)i (DATA2<N2,K) iK«iiMmi» 

;P(DaTa2(N2.1).EJ.PLAG)G:  TO  4. 

N2  =  ''  2*1 
GO  T^  3 

4  \2  =  '!2-l 
CO  TD  199 

C  **ITE  FIRST  HEADING 

199       yRlTE(6il?5)(IHE^0M)iI  =  1.16) 

105  FORMAT(lH1.16A5) 
-RlTE(6,lClft)(LABFL(I>.I=l.l2) 

106  F0RMAT(iH  ,8X. A2.4(4X. A4) ) 
103  F0RMAT(ifcA5) 

CO  112  I«1.N1 

*R  I  TEC  A, 9999)  ISA""P1<  1)«<0*TAKI»K)  iK«li*hn) 
112  CONTINUE 

C    *RlTE    SECOND    *FADPG 

wRITE(6.1?7> (U^EAOt  I  )  . I«l. 16) 
1     7  rns«iT  M  X  .  1  AAm 
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00  113  I-1.N2 

wRIT£(6.9999)I?AMP2(I)i ( DATA2 (  I  . K )  ,  K»l , MmM ) 

113      CONTINUE 

C  *EGt*  LOOP  THOUGH  E*C^  ELEMENT 

CO  6  J=1,H.MM 
C  :£T  COUNTERS  TO  ZE«0 

J 1  x  1 
C  *OVE  fLEHEMT  BEING  poOCESSED  FROM  MAIN  ARRAYS  TO  WORKING  ARRAYS, 

C  ELIMINATING  BLANKS  anO  LESS  TWAN  VALUES 

00  7  N*l.\i 

If  (OATAKNiJJ  .LE.3.C)  CO  TO  7 

XQaTH  Jl>="iATAl(Ni  J> 

J1*J1*1 

7  CONTINUE 
J1*J1~1 

J2rl 

CO  8  N=l,N2 

IF(DATA2<NiJ) .LE.3.3)  go-  tQ  a 
X0aT2(J2)=0ATa2(vj,  J)- 
*   J2*J2*1 

8  CONTINUE 
J2=J2-1 

C  CHECK  TO  SEE  WWICH  WORKING  ARRAY  IS  SMALLEST 

xchg»o.2i 

C  IF   OAtViS  SMALLEST,  INTERCHANGE  DATA  FOR  COMPUTATIONAL  CONVENIENCE 

C  LATER 

DO  10  N=i,Jl 

TEmP«X3AT21N) 

X0aT2(N)=X0ATKN) 

X0aT1(N)=TEMP 

12  CONTINUE 

TEhPsJ2 

J2sJ1 
Jl=TEMP 

XCHGsXCHG*1.3 
C  CHECK  TO  BE  SURE  SMALLEST  ARRAY  HAS  AT  LEAST  2  VALUES  OR  COGITATION 

c  .jill  Blow  up 

IF(Jl.GT.l)  GO  TO  9 
■-PlTEt6.il?)  LABEL(J) 
liC  F0RMAT(1H3,22hI\SUFFICIENT  DATA  FOR  ,A4tl3H  TO  COMPUTE  T) 

GO  TO  6 
C  CHECK  LOG  OPTION  A^Q  CONVERT  To  LOGS  IF  NECESSaRY 
9    I<"<  IQPTCJ)  ,NE.3hL0G)  GO  TO  11 
CO  12  N=i , Jl 

12  xCATlfN)=ALOGU  (XDATK*:)  > 
CO  13  N=i,J2 

13  *OAT2r4)=ALOGi:(xOAT2(N)> 

C  SET  TJkS  t0  iZ°0 
11  SUkUs?,3 

SU*U2"Z.? 

SU"X1»*.? 

SUmX2*^3.^ 

XJl*FLOAT(Jl> 

XJ?«FLOAT( J2) 

t»0OT«32RT(XJl/Xj2) 
C    CO^PlTe    SlHS    OF    TWE    OATA    anO    SUM    u 

00    14    w»i,Jl  

SUhUx3UMu*«OaT1(n)-(XOaT2CJ)»R00T) 
<;UmU2«?Umu2*(  (XDiTl(fn-(X0AT2(N)»R00T)  )*«2) 
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SUmXi=SUmX1*XDati<\> 

14  C3NTI*,JE 

30  15  \*3,J2 

SUMX2«SUMX2*X'3*T2<%) 

15  COMp.-JE 

CO'4PjTe  G 

3»(Xji«SL'MU2)-{SUHU»»?) 
COMPUTE  A'O  PRINT  ^EAN  VALUES 

XMEAlsSUMXl/XJl 

.  .  XMEA2=SU*X2/XJ2 
XVALl^XMEAl 
XVAL2sXHrA2 

IF  <XOG.ES.Z.:)  GO  TO  149 
XVAI*X^EA1 
XVAL1=XMEA2 

XVAL2«XVAL 

149  XAVEl=XVALl 
XAV£2  =  XVA|_2 

IFcrcPT(j>,NE.3HL03)  GO  TO  153 
XAVEl  =  l2.*MXAVEl) 
XAVF2«13,»»(XAVE2) 

150  i-t?lTE(6,l5l)  Li3EL{J)  .XAVEl 

151  F0RMaT{im3,i5h  "EAN  VALUE  OF  ,A4,13H  FOR  SET  i  =  ,F1Z  4) 
w"lTE(6,152)  LiBEL(J).XAVE2 

152  F0PHaT(ihz,i5m  VFA-'-J  VALUE  OF  .A4.13H  FOR  SET  2  =  .F1Z.4) 

T»(XMEA1-XVEA2J/<S3RTC3/((XJ1»«2)*(XJ1.1)))) 

T«ABS(T) 

vOEGFsJl-i 

•-«ITe(6.129)  L*3EL(J).^0EGF»T 
129  F0R»*AT(iwa,i3H  T  VALUE  FCR  ,A4,6H  kItH  »!3,22h  DEGREES  CF  FREEDOM 
'  1*  .Fi:,4) 
6  eONTiVJE 
1  C3\Tp.'JE 

E"JO 


0 
G 
0 
D 
G 
0 
0 
0 
D 
G 
D 


G 


G 


T2O09       Table  18  .  Program  used  to  Calculate  yg 

Geometric  Means  and  Deviations. 

3 1  HENS  I  OK  IFnTCift) , LABEL < 22 > , FLAG < 1* > , *LMT ( 22) , 0L< 22 ) , OATA ( 999 ,23) 
l,XnATA{999),!HEA',(26)»I0PT(22),  IVJ"(299) 
C  *Ea3  In  NU^BE*  Or  ^aTA  SEtS,  NUMBER  or  ELMENTS 
«C*0(13ill4JM»tttt«iH 

114  r-!SMAT(2!2) 

C   sJEaO    Is   Elements  OETERMnEO 

»CA0(t3il8?)(LAS£L(!)»l«i»HMM) 
122    Mrhav(12a5) 
C   REaD    In   L03   nPTin" 

5F.A3Cl3.ll7) IICPT(!),I«1#MMH) 
117    F3RMAT(12A3) 

C  REad  In  FUG  f?R  El0  CT  3ATA  SET 

REa0<13,122) jrLA5<I>iI«l*MHM) 

102  F0RMaT(12F3,:) 

C  REaO  In  DETECTION  L I " I T  FOR  EACH  ELEMENT 
REa0(13,123)(0L(I)»!s1»MmM) 

103  F0R^AT(12F5.?  ) 

ricTRsALocti:. > 

00  12  h  =  1i-*M 
C  READ  IN  HEADING  FIR  each  DATA  SET 
sEaD(13,!15) (IheaO(I) • I =1 # 16  ) 

115  FCRHAT(16A5) 
WRITE(6.H6)  (  IuEA0(  !  >  »  I«l#16) 

116  fORMA'f  (1H1.16A5) 
wRlTE(6.l2n)  <LA  =  ElU>iI  =  l,MMM) 

122     F0R»:*T(iX,i2(5»»A5)) 

N«l 
C  READ  DATA  CHECING  FOR  LAST  SAMFLE  (CAUTI0N-9999MIGHT  HAVE  TO  BE 
C  CHANGED  FOR  DATA  IN  A  DIFFERENT  FORMAT, 
5  RE»0(8i2Z:>INlT-(N)i  <DATA(\.K),KalfMMMj 

IF<INUM(V) ,£0,999)  GC  TQ  7 
2vT2  ,    F0R"AT(I3.51X.F5.1#2X,F4.2//) 
WRITE (6, 221 >  (~ATA(N4,K),K=1,9) 
2J1      F0RuAT(lX,l^(F7.2,4X)) 
C  COUNT  DATA 
N  =  N*1 
GO  T0  5 

c  begin  loop  through  each  element 

7  03  3  J=1»M«M 
C  SET  COUNTERS  Av0  SU^S  TO  HERO 
S'JfXsC. 
SUMX2»3, 
\PsZ 

••■.'■a 
•  «i 

C  *ATC-<  FOR  FlaGi  SAMPLES  BELOW  OETECTlON  LIMIT 

6  if<oata(k.,  j)  ,e:.flag<  J)  )  gj  to  1 

IF(TA*A(\,J),LT.T.  )  GO  TQ  3 
If  (  3AT A  Pi  1 J ),EQ,?»3n  GO  TO  2 

IF( lCnT( J) ,NE.3hL0G)  Gj  TO  13 
C  "3\.V£RT  TO  LnGS 

X*ATAP.»)sAt3Gi<  rOATACNi  J)) 
GO  TO  14 

13  X^aTa(n)=CATa(\, J) 

C  2lJM  nATA(LOGS)  AN3  UPOiTE 'COU\TERS 

14  SUmXsS  jmx**DaTa  (•.) 
T>X2«$U!-X2*x0ata(  ;>»»2 

G?  TC  4  ■  

3  P«MP*1 

4  T  s  H  T  *  1 


i 
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2  v«N*l 
Gn  TC  6 
C  CHEcTtO  '.£  SURE  S^mE  3^  T*E  DETERMINATIONS  ARE  ABOVl  THE  DETECTION 

C  LI-IT 

1  ir<MT.\E.D  GO  to  iP 
11  wSlTCt6»llD  LABELU)  *.-„• 

ill  rORHAT<lHZ,33H.\o  ANALYTICAL  DETERMINATIONS  FOR  , A4.21HB0VE  DETECT 
HON  LPtTJ 
GO  T0  9 
C  3ETCSJMIME  THE  VJMRER  QF  DETERMINATIONS  ABOVE  ThE  DETECTICK  LlMT 

Xi    \-JsNT-\P 

If(^iN.EQ.?J  CO  to  11 

c  co*v£rt  m'^ber  to  floating  point 

xn=fldat<nn> 
c  -rite  sums 

W*ITE<6,112>  LA3EL(JJ»SUMX 

112  FnRMAT(lue»15-^Su^(LCG)  X  FOR  .A4.3H  «  ,F12,3) 
JRITE(6,113)  LA3ELU)  .SUMX2 

113  rORMATfiH  ,l3^SuuaOG)  X»»2  TOR  ,A4,3H  =  >F12.3) 
C  DETERMINE  GEOMETRIC  ?-EaN 

XMEAP=5U^X/X\ 

C  IF  LOG  OPTION  *0T  Take**  GO  TO  BRANCH 

ir(ICPT(J),NE.3HL0G)  GO  TO  15 

GMEAPsEXPCFACT^^FAP)  -  ■ 

WRITE(6.1Z4)  L*3EL(J) .CMEAP 

104  FORMaT(ih  ,1«5-JGE^mEtR1C  MEAN  FOR  .A4.3H  •  .F13.3) 

C  DETERMINE  GEOMETRIC  nEV!ATION 

G0EV  =  EXP(FACTR»(t?C.^T((l./((XN*»2)'-XNJ)«(((XN»$uMX2)-(SufX««2)n))) 

»JRITE(6i125)  LA3EL(J)  »GDEV 
135  FORMATdH  ,24WGE0M£TRIC  DEVIATION  FOR  ,A4,3H  a  ,Fl3,3) 
C  DETERMINE  SIGMA  PRIme  (USGS  PROF.  PAPER  574-3»  P.  7) 

SP  =  SQRT(  <SUM*2/XN!>-<XMEAP»»2)> 
SPp=EXP(FACTR«SP) 
WRITE(6,106)  LA3ELU).SP 

106  FORMAT^  1I6HSIGHA  PRIME  FOR  ,A4»3H  •  .F10.3) 

C  DETERMINE  NUMBER  OF"  SAMPLES  SElO'h  THE  DETECTION  LIMIT  aNC  PROFORTION 
C  ABOVE  THE  OETE^TIO'v  L  I  m  I  T 

-RlTE£6,107)  NP,NT,LA3EL(J> 

107  F0RMAT(1H  ,13, 27^  SAMPLES  OUT  OF  A  TOTAL  OF  .13. 
135H  ARE  BELOW  The  DETECTION  LIMIT  FOR  ,A4> 

XHsFLlAT^ej/rLOaTtNT) 
»=»lTE(6il28)  LA3EL(J).XH 
138  FORMaT(im  ,46HFP.aCT!CN  OF  SAMPLES  BELOw  DETECTION  LIMT  FCR  ,A4,3H 

C  3ETE*MINE  X-aXIS  VALUE  FOR  ESTIMATING  LAMBDA  (USGS  PRCF,  PAPER  574-B, 

xlam=(sp»»?)/((x^eaR-alogi:(dl(J)))»#2) 

vRlTE(6»lir)    L»brL(J) .XLAM 
11/J    F:RMAT(1H    ,42-^ALUE    Or    X-AXIS    FOR    ESTIMATING    LAMQCA    FCR    |A4»3H    *    , 

1F1.-.3J 
GO    TC    3 
C    -JRjTr    ARlT-iMrTlc    Mf*N 

15    i-.RlTEC6.ll8)    L*BEL(  J>  .XMEA? 
lit?    FORMat<1h    ,2^-^ariThmEtIC    MEaij    FOR    .A4.3H    *    .Fi0.3> 

C  CO^PJTf  standard  deviation 

STrEv»SQPT(  (i./(  (X*<**2)*XN))«(  (  (xN«SUMX2)«(SUMX»«2)))) 

»Rl  TE(6.1l5 J    L-BEL ( J) .GTDEV  " 

119   f^omat(iw    ,23-tSTANCARH   DEVIATION   Fe*    ,a4,3h    *    ,Fl0.3]i 
SPsSCRTC  (SL^X2/X'.  )-{XMFAP»«2)  ) 
-3ITEf6.1?6)    LAfarLU).iP 
v.»|TEt6«107)    NF,\T#LA3LLU) 
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V.*ITEC6.1Z8>  LABEL(J)»XH 
C  COMPUTE  X-AXIS  VAL'uE  E3R  ESTIMATING  LAMSOA 
XLAM=(SP»*2)/(  (x^Eap-DKJ)  >»«2) 
wRITE(6»11Z)  LAHELf JJiXLAM 
9  CONTINUE 
12  CONTINUE 
END 
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Fisrure  6.   Plot  of  Analysis  of  Variance 
Sample  Locations. 
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Table  19. 


PLC 

C  "P 

j  f  ■  ■    t 

NU* 

«    M 

1 

1J1 

2 

-.  '2 

•x 

1,3 

4 

i  :4 

5 

i:6 

6 

L'7 

7 

i;s 

8 

i:9 

9 

27: 

12 

271 

11 

272 

12 

273 

13 

283 

14 

281 

15 

282 

16 

233 

17 

322 

18 

323 

19 

324 

93 

23 

325 

92 

21 

333 

9J 

22 

531 

9i 

23 

332 

9J 

24 

333 

1T9 

25 

4l6 

23 

26 

417 

98 

27 

4:8 

93 

23 

419 

9J 

29 

422 

'_-T.      LC*'.  "0 

3T5H1*.  1"5>'47J.  -?:■. 

3;5;i*.  I'S: *T3.  7.  . 

3:5-2?.  U9r3J4.  r?3, 

3-55?4.  Ild*-3?1.  4'. 

3*5r.2/.  1292313.  64. 

3»5^lP.  1252311.  59. 

3-ori:?.  1ZS7316.  51. 

2-i5b±T  .    1237316.  64. 

•355714.  1337812.  64, 

3?571?.  1332529.  93. 

395711.  133273*.  35. 

395717..  1232728.  49. 

39562?.  1282723.  31. 

395624.  1!6?73:.  <35. 

39573?.  1162733.  33. 

39573::.  1382731.  57. 

335335.  103?916."  PI. 
395*21.  1292916.  35. 
3-54Ji.  123?913.  55. 
355333.  17S3CT3.  56. 
3-5425.  135?9*>3.'  53. 
395<27.  135:  9? 1.  49. 
3 ',5*  19.  i ■'3- 9.-1.  =5. 
3  *.  3  "•  2 '.  .  1 7  8  ;■  5  :  3  .  *  2  . 
j;5*2?.  I'S;*:!''.  5?. 
3^3.^24.  139*1^9.  39. 
3v5r3'.  1  *  2  * 1 T  7  ,    *4. 

3;5*>3r.  :'8^:"4.   r7. 

3.5'7?.  1  *3  ::P9.    7.. 


Analysis  of  Variance  Data  . 

I  .  lI  J^u-C  P-t   H0SG  J  SG 

52.  l4.  .42  3.2   -ri3  27.2 

4s.  17.  1.7C  3..*   ..i6  28 .2 

*-.  17.  .f,4  7.1   .33  25 .2 

42.  17.  ,76  7.2   .61  32.2 

5'^.  15.  .74  7.7   .;j3  28.2 

6J.  15.  .46  7.3   .93  23. Z 

46.  12.  .45  7.3   .75  27. 2 
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53 


15.   .53  7.5   .76  24.2 

22.  1.62  7.9   .43  28.2 

6=.  21.  1.2*1  7.3   .96  32.2 

36.  15.  1.93  7.7   .31  3.2 

41.  12.  1.12  7.4   .74  38.7 

73.  21.  1.33  7.7   .82  28.2 

77.  21.  1.22  7.7  1.C3  33.2 

49.  17.   .33  7.4   .37  32.2 

61.  12.   .27  7.9   .52  39.2 

75.  22.   .63  7.3  1»23  33.7 

66.  21.  1.22  7.4  l.l'J  37.2 

51.  24.  1.23  7.2  1.23  59.2 

72,  19.  1.42  3.1   .46  29.1 

6.:.  '7.  1.42  9.3  2-33  35.  C 

7...  47.   ,<$1  3.3  1.33  34.2 

•i '! .  16.  3 .  :  I  7.9   .91  33.2 

4  .  17.  2.7J  7.4   .£3  32.3 

6J.  l».  1.72  7.7   .33  2S.2 

57.  22.  1.72  7.5   oJ  26.2 

6c.  31.  l.hJ  7 . § 


^  ■   .  » 


•  7 ■   31.2 


A"*. 


2.1. 


7.4    .95   J/.2 


83 

e 

161. ? 

163.2 

132.3 
123.3 

139.3 

143.2 

126.:' 

98.9 
125.3 
221.3 

89.3 

86.1 
169.3 
293.3 
137.3 

63.7 
132.3 
152.2 
132.3 
169.3 
112.2 
145.2 
141.3 

56.2 
176.3 
172.3 
126.  J 
147.2 
14  3.3 


1.8u 

1.4  J 
,82 
1.9J 
1.8  J 
1.13 
2,23 
1.43 
1.1  J 
2.13 
1.22 
1.53 
3. 23 
3.03 
1.53 
2. 7  j 
2.23 
2,13 

2.8;' 

1.5  J 
4.7:: 
5.2.) 
,7.j 
1.5i 
1.4.1 
1,53 
2.8.1 
4.8.) 
1.9  J 
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32      4?3        3*5^,24.    Ua?l32.         "1.      7i.      27.    1,33      7.1      .553      39.2      163.3      2.23 

31  424         3?5&2*.    13S?2:6.         48.       57.       16.    2.62      7.5    1.4J      37. Z       133. J      2.33 

32  *?5        395&2i.    1383229.         96.       ft*..       1?.    2.e3      7.6    1-2'J      34. Z         88.9      1.93 


0 

1 
ffl 
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Table  20.   Analysis  of  Variance  Data 
for  Components  not  Analyzed  on  the 
Grid  Samples. 


Samp 

SG 

SG 

W, 

3G 

WG 

4(i 

WG 

R2- 

*G 

SOIL 

MUM 

Z\ 

H>*. 

Mr, 

ZN 

ZN 

8 

MS 

H 

MQ 

AS 

121 

3 

29 

2 

11 

.75 

6 

1*2 

2 

11 

.39 

6 

123 

2 

8 

.83 

6 

l'w'4 

3 

32 

2 

12 

.9j 

7 

126 

? 

4.; 

2 

9 

1.5 

8 

127 

1 

31 

12 

8 

.69 

6 

128 

5 

33 

9 

11 

.68 

6 

109 

4 

7 

.73 

4 

270 

5 

-2* 

34 

7 

9 

5 

.89 

12 

1.3 

? 

271 

3 

23 

23 

2 

11 

17 

1.1 

13 

1.2 

9 

272 

5 

22 

1 

13 

8 

?73 

4 

7 

25 

.47 

14 

1.1 

ih 

282 

3 

28 

47 

1 

7 

15 

1.1 

9 

1.5 

6 

281 

3 

34 

29 

2 

12 

24 

1.7 

16 

1.1 

8 

282 

2 

-22 

127 

11 

9 

15 

.33 

11 

1.1 

3 

283 

8 

34 

23 

.4 

11 

13 

1.3 

8 

.83 

7 

323 

5 

44 

47 

2 

9 

11 

1.4 

13 

1.4 

lfl 

324 

4 

34 

29 

1 

8 

13 

1.4 

9 

1.4 

10 

326 

3 

-2  c1 

21 

11 

7 

17 

.76 

27 

.93 

8 

327 

2 

24 

29 

4 

9 

16 

.75 

16 

1.2 

8 

330 

7 

-22 

39 

3 

14 

18 

1.6 

9 

2.0 

16 

331 

.3 

21 

-23 

4 

11 

m 

1.4 

11 

1.4 

14 

332 

2 

46 

-2" 

1 

13 

15 

.68 

11 

1.2 

6 

333 

2 

-22 

52 

11 

6 

12 

1.2 

12 

7" 

7 

416 

.4 

-2: 

-22 

12 

15 

13 

15 

.83 

12 

417 

4 

21 

-2  ' 

1 

11 

15 

.83 

11 

.79 

I? 

418 

s 

-2  ? 

-2: 

2 

7 

16 

1.3 

6.7 

1.5 

11 

419 

3 

-22 

-2  2 

6 

14 

15 

1.2 

6.7 

1.5 

15 

422 

3 

25 

5  J 

2 

8 

22 

1.7 

9.3 

1.4 

12 

423 

4 
J. 

24 

4t 

1 

12 

17 

1.5 

9.3 

.77 

13 

424 

2 

-2  J 

9 

4 

17 

1.5 

9.3 

1.1 

11 

425 

.3 

43 

-23 

2 

11 

14 

1.1 

3 

.7-J 

12 

I 

D 

0 
D 
0 
0 
0 
0 
0 

e 

0 
0 
D 
D 
I 
I 
1 
0 
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I 
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PLOT  OF  TRACT   C-A  GRID  SAMPLE: 

?19    223    221     212    236    227    233    118    117    114    37      36      31       33      23      19      21 
-„»«»■*■■■■■■■■■■ 

J  '13    2lS    214    211     239    2:3    123    119    41       40      34      33      23      29      16      17      18 

?i"    »»»■■»»■■*■■■■■■ 

1217    213    163    167    162    163    1*1     122    43      42     33     32      27      26      15      14      13 

J: 

*1   173    171     169    166    161     164    133    137    56      55     39     33      24      2S      10      11      12 

?J  '16    172    173    165    159    160    139    136    57      54      52      S3      23      22      2        3        4 
ftlri        «■■*■■■■■■■■■■■ 

•   131    177    176    175    1S3    156    140    135    61       S3      51      S3      46      *7      1         9        5 

g<i      •■•«»      »■■■■■      ■■•■■ 

ft- 

:    130    179    173    17*    157    143    141     H2    60      59      45      44      49      *8      3        7        6 

§»  M  »■»»■»■■*  ■■■  " 

132  133  136  18*  155  15*  133  13*  113  115  103  90   37   9*   65   62  63 


ll         ■         ■         ■ 


&: 


133    195    191     192    153    152    131     132    111    116    181    89      36      83      64      67      66 
I  137    190    19*    193    151     153    130    129    112    110    132   91      99      85      70      69      63 


S"   13?    197    196    195    143    1*9    127    123    139    137    96      95     93  _92     71      72     73 
;,        ■•«■!■*■        ■■■■■*■•■ 

J: _    ,y,    ,,.     iya    ,4A    143    yZb    ,33    133    136    133    97      5a      74      73      79      83 

&h      » "    « "    «     "■"■»«■«■■■■■■■ 


i;    i-i    ITS    135    1»1    ??      93      81       37      75      76      77 


■  .  M  r - -*t »- ■  ■      T~»--r1 r- p« ■ 


T*- 


IM.JJ 


LONGITUDE 


Figure  5.   Plot  of  Grid  Sample  Locations. 


*2tKW  Table  21.   Data  from  Grid  Samples.  °' 

as: 5^: LiT-    LCN-   "a  s.  &  t  PH  H°"  >"  m«  & 

i  :c:  395511.  1237415.  47.  a:.  21.  .69  8.2  .63  35.6  148.2  t.12 

2  111  395526.  1262415.   74.   57.  14.   .50  7.9  .53  31.3  124.fi  1.3Z 

3  112  395526.  1237331.  ;27.   56.  13.   .65  7.9  .64  31. 6  113.2  1.52 

4  n3  395526.  1337311.   37.   5-.  14.   .53  7.9  .82  3.' .  5  119.2  1-52 

5  114  395511.  1282311.   49.  263.  13.   .46  7.9  .83  iS.7  122.2  1.52 

6  115-  355431.  1332311.   39.   52.  13.1.13  7.6  .42  19.4  97.1  l.i8 

7  n6  395,31.    1232331.    .45.       84.  2i.    1.33  8.2  .85  25.3  126.2  1.72 

8  117  395431.    1367415.    -2?.       62.  15.    1.32  8.'2  .6J  29.2  US. 2  1.32 

9  113  375511.1232331.       42.       81.  16.       .65  8,2       .23  36.1  151.2  1.92 
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IIS  119  395626.    1282415.'    -2'..  73.  16.    1 .  IS  8.2  .38  26.7  176.2  Li2 

11  122  395626.    1782331.'    -2?.  71.  17.       .66  8.1  .84  26. J  H5.2  1-42 

12  121  395626.  1232311.   85.  62.  15'.   .65  S.'l  .«6  35.7  126.2  1.52 

13  123  395622.  1282311.  34.  62.  14.  .65  7.9  .73  33.6  127.2  1.82 
i4  124  395622.  1282331V  79.  66.  15,  1.23  7.9  .46  33.6  141.2  1.08 
J  126  395622.  1362415.   22.  69.  19.  1.92  8.3  .87  32.1  148.2  .83 

16  127  395722.  1332415.   36.  66.  17.   .73  8.3  .79  31.7  162.2  2.22 

17  128  395722.  1232331.   23.  112.  13.  1.27  SV2  .47  29„5  118.2  1.42 

18  129  .  39572?.  1332311.  24.  61.  16.1.13  8,1  .67  23.7  119.2  1.52 
It  134  395717.  1232331.  25.  79.  24.  2.42  8Vl  .55  29.9  125.2  1.62 
23  135  3^5717.  1182415.'   32.  223.  16.   .69  8.'2  ,8J,  34.3  .88.2  1.52 


21   136 


3^5717.  1232311.   59.   82.   19.  1.32   8.'3   .33   42. 9   127.2    .75 

22  139    395526.  1332435.'  -2*'.   52.   13.   .77   S'.'l   .66   27.5    56.1   2.22 

23  142    375S26.  1182519.  -2*.   69.   22.  1.12   8.3   .44   22.7   152. 2   1.22 


i 


24  143  3956'*..  1292519.   24.  114.  16.   .88  3.1  1.32  33.7  58.  C  .32 

25  :<4  3156-6.  .153?««.  -2f.  33.  14,   .25  3.3  1,5j  35.'.  123.2  1.52 

26  145  39562?.  1332435.  -2'.  72.  17.   ,93  3.1  ,39  23.4  117.2  Ll2 

27  146  39562?.  1132519.   35.  72.  13.  1.32  8.2  .46  22.4  HZ. 2  1.22 

28  147  3357*?.  1297519.  -2%  64.  17.   .57  t'.'l  .63  1?.?  144. C  1.88 

29  148  3^57-?.  1255435*;   24.  7T.  2.'.   ,3l  7.4  .*2  28,3  41.3  Ll2 
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G 
D 
0 
0 
0 
Q 
0 
0 
I 
0 
1 
1 
I 
I 
i 
1 
I 
I 
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335717.  1?S?*35,  -2".  62.  14.   .35  8.1  .67  28.8  121, 2  .32 

395717.  1?8?519."   47.  31.  23.  1.3r  8.'2  .77  29.3  163.2  1.52 

3-J5a27.  138~623.   S«.  64.  lb.   .88  8,1  ,67  23.2  122.  2  2. 22 

335*22.  1132623.   47.  94.  15.   .65  8,1  .59  ?5 . 4  47.4  .92 

3v573?.  1382623.   33.  133.  22.  1.73  7.9  .37  22.7  59.2  1.22 

3937??.  13326Z3.   23.  93.  23.   .66  8.2  .83  32.4  146.2  L32 

395717.  133?6"3.   34.  8-".  25.  1.3a  7.9  ,65  32.0  125.2  1.42 

395717.  138P623.   29,  53.  16.   .89  8,2  ,3l  27.5  54.4  1,12 

3956J6.  13826:3.'  -2?.  35.  14.   .68  8.3  ,9a  34,1  124.2  .82 

^395626.  12S2623.  -23.  59,  15.   .86  8.'2  l.U  33.5  53.5  1.22 

39572?.  1332728.   24.  75.  23.   .75  8.2  .73  38.2  161.2  2-l2 

39573?.  1782723.  -23.  66.  17.  1.23  7.'9  .63  3C .  5  gfi.i  ltle 

39562?.  12827?3.  -2?,  31.  15.  1.2.1  8,3  ,4i  34,1  m.2  i,72 

395622.  1352723.  -23.  235,  2l.  1.2.T  8'.2  .45  34.5  143.2  1.22 

395431.  12826?3.   48,  46.  16.  1.53  7','9  .55  27.3  53.4  .52 

39543j.  1382623.   35.  72.  19.  1.2a  3.'2  ,57  28.7  123.2  1.10 

395511.  1332519.   83.  72.  15.  1.9J  8*,0  ,47  24.4  128.2  1.22 

395511.  1282435.   53.  78,  -17.  1.6C  7.9  ,42  26. J  121.2  1-82 

395431.  1332435.   23,  66.  13.   .94  3.2  .25  32.5  114.2  1,02 

395431.  1232519.'   3?.  68.  17.   .77  8.1  1,32  23.5  155.2  1.12 

395511.  12826?3.  129.  4S,  13.  1.2C  8.2  .36  36. 1  123.2  1.82 

395511.  1337623.   6A.  57.  13.  1.62  8.1  .43  39.; 3  55,3  lt2e 

355526.  133P623.   3?.'  65.  16.  1.53  8.2  .43  3'« .  1  £2.2  .32 

39'526.  1*S?6?3.   52.  55.  17.   ,76  8,2  ,3i  ?6.!3  122.2  1.22 

3-;5526.  1382778.  -27.  7*6.  3J.  3.3.:  8.2  .23  2. J  123.2  2.32 

39563*.  Ii3"77d.  -2'.  53.  Id.  1.22  7,8  .49  24. d  127.2  1,32 

395^.^.  1282728'.  -2*.  214.  19.  1.7.1  8.1  ,5j  26.6  117.2  1.82 

395526.  1232723.   2ft.  73.  19.   .83  8',**  .4j  3C.5  1«6.2  1,88 

395511.  I.!3?7e8.   29.  7-.  14.  1.5,'  8.1  .61  36. 3  119.2  1.32 

3^5*31.  l*2?7ra.   3'J.  ?5.  2-.  1.2:'  7.3  .63  3^.3  132.2  1.12 
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?39      4.8         3^57'?.    l,83C'-3.       2?.       70.       1^.     i.i.«      7.7       .3d      29.6       124.2       1-60 
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8 
212   4H    3957.1?.  1287916.   37.  162.   15.   .63   8,3   ,55   32.5   124.2    .72 
8 

211  4i2    39572?.  1283222.'  -22.   66.   15.   .86   7,9   ,4i   34.1   122.2   Lie 
92 

212  4i3    395717.  1283222.   56.   89.   14.  1.10   7.5   ,5j   37.6    48.2   l.lZ 
8 

213  414    39562?.  1287124.'   91.   77.   37.  2.22   8.2  1.3j   24.9   167.2   2.32 
98 

214  415    39573?.  1232124."   6?.   72.   32.  3.22   8,1   .64   36.3   128.2   2.2Z 
93 

215  422    39572?.  1252:24.   5.7 .   55.   3d.  2.3C   8,*1   .74   42.4   124.2   2.32 
3 

216  421    395526.  1263279.   2*.   72.   19.   .75   8.2   .78   36,7   122.2   2.62 
92 

217  4?6    39562?.  1",S3229.   25.   66.   24.  2.32   7,8   .45   24.8   114.2   2.82 
92 

218  427    39572?.  1232229.   26.   61.   34.  2.22   8.1   .82   35.9   128. 2   4»i2 

0 

219  428    395717.  1283229.   37.   72.   14.  1.6fl   7,9   ,43   28.9   126.2   1.52 
90 

222   429    395717.  1283124.   24.   65,   59.  1.22   8.2  1.43   26,8   152.2   1.22 

92 
221   433    395717.  1?83L?4.   62.   61.   56,  2.52   8.1  1.52   33.3   112.2   3.92 


T20O9 


95 


PLOT  OF   GROUP  i  UINTA  FORM. 
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Figure  18.   Flot  of  Group  1  Uinta  Formation 
Samole  Locations. 
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Figure   19.      Plot   of  Group  2   Ulnt-5    Formation 
SairtJle    Locations, 


T2009 


97 


PLOT  OF  GROUP  3   UINTA   FORM. 
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Figure  20.   Plot  of  Group  3  Uinta  Formation 
Sample  Locations. 
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